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Dynamic  Multi-axial  Loading  Response  and  Constitutive/Damage  Modeling  of 

Titanium  and  Titanium  Alloys 


Principal  Investigator:  Akhtar  S.  Khan 
Department  of  Mechanical  Engineering 
University  of  Maryland,  Baltimore  County 


Research  Objectives 

Titanium  alloys,  and  ceramics  encapsulated  in  titanium  alloys  are  considered  for  next 
generation  of  lightweight  tanks.  The  objectives  of  this  research  project  are  to  obtain 
responses  of  newly  developed  economical,  but  with  much  higher  Oxygen  content,  Ti- 
6A1-4V  alloys,  determine  the  suitability  of  existing  constitutive  relations  to  model  these 
responses,  and  to  detennine  material  constants  for  these  models  (e.g.,  JC  and  KHL 
models)  for  incorporation  in  computer  codes  for  penetration  simulation. 


These  objectives  are  being  achieved  through  uniaxial  loading  experiments  over  a  wide 
range  of  strain  rates  and  temperatures  in  the  first  year  of  funding,  and  through  multi-axial 
loading  experiments  over  a  wide  strain  rate  range  and  a  wide  temperature  range  in  the 
second  and  third  year  of  funding,  respectively. 

Approach 

Uniaxial  compressive  experiments,  over  a  wide  range  of  strain  rates  and  temperatures,  are 
used  to  determine  material  constants  for  two  constitutive  models  (Johnson-Cook  and 
Khan-Huang-Liang);  experiments  are  perfonned  on  several  titanium  (Ti-6A1-4V)  alloys. 
The  results  from  these  models  are  correlated  to  these  experimental  data  to  demonstrate 
the  flexibility  in  each  model.  Most  other  investigators  who  have  performed  similar 
studies,  stopped  at  this  stage.  The  approach  taken  by  other  scientists  is  not  right  because 
if  correlations  are  done  with  the  same  experiments  from  which  the  material  constants  are 
determined,  then  correlations  should  be  always  good  for  any  constitutive  model.  We  have 
gone  one  step  further,  which  is  very  desirable,  if  not  required,  to  validate  a  constitutive 
model.  In  the  phase  of  the  current  year  of  funding,  multi-axial  loading  experiments  (e.g. 
dynamic  torsion  followed  by  dynamic  compression  using  Kolsky  or  split  Hopkinson  bar 
technique,  uniaxial  compression  followed  by  biaxial  compression,  etc)  have  been 


performed  and  compared  to  predictions  from  these  two  constitutive  models  to  establish 
validity  of  these  two  models  over  a  wide  ranges  of  strain  rates  and  temperatures. 


Background 

Since  the  introduction  of  titanium  and  titanium  alloys  around  1950,  they  have 
become  important  materials  for  aerospace,  energy,  and  chemical  industries.  They  are 
used  chiefly  for  parts  that  require  good  corrosion  resistance,  moderate  strength  up  to  588 
K,  and  lightweight.  Titanium  alloys,  especially  Ti-6A1-4V,  an  a+P  type  titanium  alloy  is 
largely  used  alloy  in  many  industries  because  of  its  extremely  attractive  properties  like 
high  specific  strength,  good  deformability,  reasonable  ductility  and  ability  to  withstand 
high  temperatures  and  resistance  to  corrosion.  It  is  primarily  used  in  aero-engine,  gas 
turbines  and  other  applications.  The  development  of  a  relatively  economical  Ti-6A1-4V 
alloy,  with  a  low  interstitial  content  has  prompted  a  lot  of  interest  in  its  possible  use  in 
annor  tanks  because  of  improved  ductility,  whereas  the  conventional  more  expensive  Ti- 
6A1-4V  alloy  has  been  used  primarily  in  aerospace  components  chiefly  because  of  its 
high  strength  to  weight  ratio.  This  advantage  has  also  been  contemplated  in  various 
applications  like  armor,  including  ceramic  tiles  encapsulated  in  these  titanium  alloys. 

However,  applications  of  this  material  especially  in  defense  structures  require  a 
thorough  understanding  of  the  mechanical  properties  and  deformation  responses  under 
different  loading  conditions.  Of  all  the  kinds  of  deformations,  dynamic  deformation  is  of 
prime  importance  when  it  comes  to  use  for  military  applications  and  it  should  be  studied 
in  detail.  Also,  the  dynamic  deformation  investigation  need  to  be  multiaxial  in  nature  so 
as  to  simulate  the  real  life  applications.  Although  multiaxial  dynamic  deformation  is  most 
essential,  one  needs  a  complete  set  of  experimental  observations  to  make  essential 
judgments  regarding  the  applicability.  Also,  one  needs  to  model  the  constitutive  response 
for  dynamic  rate  of  loadings  as  well.  The  constitutive  model  which  is  required  for  finite 
element  and  hydro-code  simulations  has  to  be  simple,  with  lesser  number  of  material 
constants  to  reduce  calculation  times  and  should  also  depend  on  experiments  which  can 
be  perfonned  with  ease  to  reduce  experimentation  time  and  thus  be  cost  effective.  The 
present  investigation  emanates  from  these  requirements. 

Different  constitutive  models  have  been  proposed  over  the  last  few  decades 
describing  the  material  behavior  for  high  strain  rate  applications  in  the  plastic  regime. 
These  can  be  classified  into  two  categories;  the  purely  phenomenological  ones,  e.g., 
Johnson-Cook  (JC)  and  Khan-Huang-Liang  (KHL)  models,  and  the  physically  based 
models  e.g.,  the  ones  by  Zerilli  and  Armstrong,  Mechanical  Threshold  Stress  and  the  one 
proposed  by  Nemat  Nasser  et  al.  The  phenomenological  models  usually  are  very  simple 
in  nature  and  have  lesser  number  of  material  constants  involved  than  the  so  called 
physically  based  models.  Also,  the  experiments  required  for  the  determination  of  the 
material  constants  in  the  phenomenological  models  are  easily  performed  than  those 
required  by  the  physically  based  models.  The  two  phenomenological  models  are 


compared  and  their  capabilities  to  model  complex  material  behavior  are  understood  in  the 
first  phase  of  this  investigation. 


Also,  it  is  a  well  known  fact  that  the  mechanical  properties  of  this  titanium  alloy 
are  heavily  influenced  by  its  impurities  and  heat  treatment,  so  it  is  imperative  that  this 
alloy  be  understood  in  terms  of  the  microstructure,  constituents,  and  the  processing 
history.  The  unalloyed  titanium  exists  in  mainly  two  forms,  a-titanium  (hep  structure)  at 
room  temperature,  and  at  high  temperatures  it  exists  as  P-titanium  (bcc  structure). 
Alloying  of  this  material  is  performed  by  adding  aluminum  and  vanadium,  whereas 
oxygen,  nitrogen,  and  carbon  are  the  interstitial  impurities.  Aluminum,  oxygen,  nitrogen, 
and  carbon  are  the  a  stabilizers,  while  vanadium,  iron  and  manganese  are  the  P 
stabilizers.  The  a  phase  (primary  and  secondary)  is  more  dominant  in  the  temperature 
ranging  from  300  K  to  800  K,  while  the  amount  of  P  phase  starts  to  increase  thereon  and 
is  almost  totally  P-phase  at  about  1270  K,  which  is  called  the  P  transus  temperature.  The 
a  phase  is  much  harder  than  the  P  phase,  so  the  drop  in  stress  levels  upon  increase  of 
temperatures  can  be  noticed  with  increasing  P  content.  However,  formability  increases  as 
the  phase  changes  from  a  to  p.  The  uniaxial  behavior  of  three  different  alloys  with 
different  alloying  and  interstitial  content  are  experimentally  studied,  modeled  and 
understood  in  terms  of  the  chemical  composition,  microstructure,  etc.,  in  the  second 
phase  of  this  investigation. 

In  the  last  phase  of  this  investigation,  the  multiaxial  response  of  one  of  the  alloys 
which  was  manufactured  using  single  electron  beam  melting  process  was  studied  and 
modeled.  Dynamic  torsion,  dynamic  torsion  followed  by  dynamic  compression 
experiments,  and  non-proportional  uni-axial  to  biaxial  experiments  were  perfonned  to 
understand  the  multi-axial  response  of  the  alloy.  Because  of  the  HCP  crystal  structure  of 
the  alloy  at  low  to  moderate  temperatures  and  production  in  a  plate  geometry,  this  alloy 
has  been  known  to  exhibit  anisotropy.  Also,  there  is  a  significant  difference  in  the  flow 
stress  in  tension  and  compression.  An  anisotropic  model  was  chosen  with  capability  to 
model  the  anisotropic  response  and  also  the  strength  differential  along  with  the  proposed 
constitutive  model  to  simulate  the  response  under  complex  dynamic  multiaxial  loadings. 


Executive  Summary  of  Achievements  I 

Dynamic  deformation  has  been  of  interest  not  only  in  impact  and  penetration 
related  problems  but  also  in  high  speed  machining.  Titanium  alloys  have  been  studied  by 
several  investigators  because  of  their  use  in  aero-engine,  gas  turbines  and  other 
applications  due  to  their  high  strength  to  weight  ratio,  ductility,  and  ability  to  withstand 
high  temperatures  and  resist  corrosion.  The  development  of  relatively  economical  Ti- 
6A1-4V  alloy,  with  resulting  high  oxygen  content,  has  sparked  interest  in  its  possible  use 
in  lightweight  tanks;  the  conventional,  more  expensive  Ti-6A1-4V  alloy  has  been  used 
primarily  in  aerospace  components.  The  potential  applications  in  armor,  including 
ceramic  tiles  encapsulated  in  titanium  alloys,  have  motivated  several  studies  (Gray  III, 


1997;  Follansbee  and  Gray,  1989;  Lesuer,  2000;  Nemat-Nasser  et  al.,  2001;  Majorell  et 
al.,  2002). 


However,  most  of  these  studies  have  been  insufficient  and  not  comprehensive 
(Follansbee  and  Gray,  1989)  to  ascertain  the  material  response  under  different  loading 
conditions,  or  the  range  of  strain  rates  in  the  investigations  has  been  insufficient  (Nemat- 
Nasser,  2001).  Also,  some  inadequacies  and  contradictions  are  found  in  the  published 
modeling  efforts.  The  constitutive  models  used  for  high  strain  rate  applications  can  be 
classified  in  two  categories;  the  purely  phenomenological  ones,  e.g.,  Johnson-Cook  (JC) 
(Johnson  and  Cook,  1983)  and  Khan-Huang-Liang  (KHL)  models  (Khan  and  Huang, 
1992;  Khan  and  Liang,  1999;  Khan  and  Zhang,  2000,  2001)  and  so-called,  “physically 
based  models”  e.g.,  the  ones  by  Zerilli  and  Armstrong  (1987),  Mecking  and  Kocks 
(1981),  etc.,  that  were  used  frequently  by  Follansbee  and  Gray  (1989),  and  in  a  modified 
form  by  Nemat-Nasser  et  al.  (2001).  The  latter  group  discusses  the  mechanisms  of  plastic 
deformation,  mainly  dislocations.  However,  the  material  constants  are  determined  by  not 
measuring  any  deformation  mechanism  related  quantity,  but  by  choosing  constants  to 
“fit”  the  uniaxial  stress-strain  curves  at  different  strain  rates  and  temperatures,  just  like 
the  purely  phenomenological  models.  Mecking  and  Kocks  model,  as  used  by  Follansbee 
and  Gray  (1989)  has  23  constants,  while  these  constants  range  from  12  to  8  in  case  of 
Cheng  and  Nemat-Nasser  (2000)  and  Nemat-Nasser  et  al.  (2001),  respectively, 
depending  on  whether  they  include  modeling  of  dynamic  strain  aging  or  not.  Johnson- 
Cook  and  Khan-Huang-Liang  models  have  5  and  6  constants,  respectively.  Johnson- 
Cook  and  Khan-Huang-Liang  models  are  used  in  this  investigation  due  to  their 
advantage  of  fewer  constants  and  their  ability  to  model  the  observed  material  response  as 
closely  as  with  models  with  many  more  constants. 


Modeling  approach 

The  material  constants  for  both,  the  modified  KHL  and  JC  models  were 
determined.  Using  these  material  constants,  correlations  were  obtained  and  compared  to 
experimental  results.  The  modified  KHL  model  is  as  follows: 
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where,  a  is  the  true  (Cauchy)  stress  and  ep  is  the  true  plastic  strain.  Tm ,  T ,  Tref  are 
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Fig.  1.  Comparisons  between  KHL  model  and  JC  model  correlations  for  alloy  3  studied  in  the  investigation 

at  room  temperature  and  different  strain  rates. 


upper  bound  strain  rate)  and  £=i  s" 1  ( re fcrcncc  strain  rate,  at  a  reference  temperature  of 
Tref,  usually  room  temperature,  at  which  material  constants  A,  B  and  no  are  determined). 
e  is  the  current  strain  rate,  n , ,  C  and  m  are  additional  material  constants.  For  Ti-6A1-4V 
alloys,  the  melting  temperature  was  taken  to  be  1933  K  (ASM  handbook,  1994).  The 
reference  temperature  was  the  constant  room  temperature  for  experiments  at  296K. 

The  JC  model  which  has  been  used  previously  is  given  as  follows: 
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The  KHL  model  has  two  distinct  advantages  over  the  JC  model.  These  are  the 
addition  of  another  material  constant,  n] ,  which  is  able  to  simulate  the  decreasing  work 
hardening  behavior  of  certain  materials  with  increase  in  strain  rate.  The  second  advantage 
is  the  modified  temperature  term  which  enables  to  simulate  the  material  behavior  below 
the  reference  temperature  Tref.  In  the  JC  model  the  temperature  term  cannot  accommodate 
the  case  when  the  current  temperature  is  lower  than  the  reference  temperature  as  then  the 
above  term  within  parenthesis  becomes  a  negative  number  raised  to  the  power  m.  The 
model  material  constants  are  determined  through  systematic  procedure  and  further 
refined  using  least  square  optimization  technique.  An  illustrative  example  of  the 
modeling  using  JC  and  KHL  is  shown  in  Fig.  1  for  one  of  the  alloys  studied  during  the 


current  investigation.  The  figure  clearly  shows  the  capability  of  the  KHL  model  against 
the  JC  model  in  predicting  the  material  behavior  under  different  loading  conditions.  Also, 
the  JC  model  and  KHL  model  comparisons  are  shown  for  a  similar  alloy  data  published 
in  the  literature  in  Fig.  2. 


Fig.  2.  KHL  model  and  JC  model  correlations  with  experimental  observations  for  experiments  performed  at 
room  temperature  [Data  from  Macdougall  and  Harding  (1999)]. 

To  determine  the  temperature  dependence  in  the  KHL  and  JC  models, 
experiments  were  performed  on  a  titanium  alloy  at  different  temperatures.  The 
temperatures  chosen  for  these  experiments  were  both  above  and  below  the  reference 
temperature  which  in  this  case  was  room  temperature,  but  below  0.4Tm.  The  experimental 
observations  along  with  the  model  correlations  are  shown  in  Fig  3.  Note  the  absence  of 
any  JC  correlation  for  the  experiment  perfonned  below  the  reference  temperature  of 
233K;  the  KHL  model  is  able  to  correlate  the  material  response  for  that  experiment  also 
quite  well.  Overall,  the  KHL  model  was  able  to  simulate  the  material  behavior  at  these 
temperatures  reasonably  well  whereas  JC  was  unable  to  simulate  for  certain 
temperatures;  also  the  correlations  obtained  through  JC  model  were  also  not  as  good  as 
the  ones  obtained  by  using  KHL. 


Fig.  3.  KHL  and  JC  model  correlations  for  experiments  performed  at  different  temperatures  for  alloy  3  at 
the  strain  rate  of  10-3  s-1.  Note  the  absence  of  correlation  at  233Kwith  JC  model  in  which  the  temperature 

term  (T-Tn,f)l(Tm-Tref)  becomes  negative. 

From  the  current  phase  of  investigation  it  is  found  that  Ti-6A1-4V  is  non-linearly 
dependent  on  strain  rate,  as  well  as  temperature.  The  KHL  model  is  found  to  correlate 
better  than  JC  model  especially  during  dynamic  deformation  regime.  The  thermal 
softening  at  high  strain  rates,  together  with  reduction  in  the  work  hardening  rate  with 
increase  in  strain  rate  and  strain,  is  captured  much  better  by  the  KHL  model  than  the  JC 
model.  The  temperature  term  of  the  modified  KHL  viscoplastic  constitutive  model  is  able 
to  correlate  well  with  the  response  at  a  temperature  (233  K)  which  is  lower  than  the 
reference  temperature  (296  K);  the  JC  model  is  not  valid  in  this  case. 


Executive  Summary  of  Achievements  II 

Results  from  a  series  of  experiments  on  three  different  titanium  alloys,  under 
quasi-static  and  dynamic  loading  conditions  are  presented.  These  alloys  have  been 
designated  with  different  numbers  i.e.  as  alloys  1,  2  and  3.  Alloys  1  &  3  were 
manufactured  using  single  electron  beam  melting  process,  resulting  in  higher  oxygen 
content  than  the  ELI  version  (alloy  2). 


The  strain-rates  varied  from  10"6  s'1  to  3378  s’1,  while  observations  are  made  at 
temperatures  from  233K  to  755K.  Differences  in  these  alloys  are  observed  in  terms  of 
thennal  softening,  work  hardening,  and  strain-rate  &  temperature  sensitivities.  The  Khan- 
Huang-Liang  (KHL)  model  is  used  to  effectively  simulate  the  observed  response 
obtained  from  these  experiments.  The  model,  with  constants  determined  from  the  above 
experiments,  is  then  used  to  predict  strain-rate  jump  experimental  results  and  also  high 
temperature  dynamic  experiments  for  one  of  the  alloys;  the  predictions  are  found  to  be 
very  close  to  the  observations.  The  chemical  composition  of  the  three  alloys  used  in  this 
investigation  is  shown  below. 


Material 

Al 

V 

Fe 

Y 

H 

N 

o 

C 

Ti 

Oeq 

Alloy  1 

6.26 

4.16 

0.14 

<0.0003 

0.0031 

0.008 

0.178 

0.047 

REM 

0.229 

Alloy  2 

6.30 

3.86 

0.18 

<0.0003 

0.0026 

0.003 

0.112 

0.045 

REM 

0.152 

Alloy  3 

5.97 

4.09 

0.15 

<0.0003 

0.0041 

0.008 

0.174 

0.043 

REM 

0.222 

Table  1.  Chemical  composition  of  the  three  Ti-6A1-4V  alloys  used 

The  equivalent  oxygen  content  (Oeq  =  O+2N+0.75C  ),  proposed  Conrad  et  al. 

(1975)  is  also  shown  in  the  table  above  which  gives  the  effect  of  dislocations-impurity 
interaction  on  the  yield  strength  of  the  material.  It  can  be  seen  clearly  that  alloy  1  and  3 
are  the  higher  oxygen  content  titanium  alloys,  while  alloy  2  is  the  ELI  version  of  the  Ti- 
6A1-4V.  The  microstructure  of  these  alloys  was  also  detennined  before  and  after 
deformation  using  a  scanning  electron  microscope  (SEM).  The  specimens  were  cut  along 
the  compression  axis  of  the  specimen  (the  thickness  direction  of  the  plate).  The  initial 
microstructure  of  one  of  the  higher  oxygen  content  alloy  (alloy  3)  and  the  ELI  version 
(alloy  2)  are  shown  in  Fig  4.  It  clearly  shows  distinct  features  in  the  microstructure  of  the 
two  Ti-6A1-4V  titanium  alloys.  Further  studies  by  Conrad  et  al  (1975)  noted  that  the 
strengthening  of  the  alloy  due  to  these  interstitials  was  dependent  on,  in  decreasing  order, 
Carbon,  Oxygen  and  followed  by  Nitrogen.  Effect  of  Hydrogen  in  most  cases  was 
neglected.  They  also  found  that  between  temperatures  from  300  K  to  800  K,  the  flow 
stress  vs.  temperature  curve  of  a-titanium  was  almost  parallel  to  that  of  Ti-6A1-4V  alloy 
of  nearly  same  interstitial  content,  concluding  that  alloying  elements  had  insignificant 
effect  on  the  temperature  sensitivity  in  this  range;  the  interaction  of  dislocations  with 
interstitial  impurities  (C,  N,  O  and  H)  had  more  effect  on  the  response. 


a) 


b) 


Fig.  4  Initial  microstructures  of  the  (a)  ELI  version  of  the  titanium  alloy  (alloy  2),  and  (b)  the  higher 
oxygen  content  titanium  alloy  made  from  single  electron  beam  melting  process  (alloy  3). 


In  the  current  investigation,  it  is  observed  from  the  measured  responses  that  alloys  1  and 
3  have  higher  flow  stress  as  compared  to  the  alloy  2.  This  is  attributed  to  the  fact  that  the 
interstitial  solute  content  or  the  equivalent  oxygen  content  ( Oeq ),  which  is  known  to 
strengthen  the  alloy,  is  higher  in  these  alloys. 


The  material  constants  for  this  model  are  determined  using  several  uniaxial  quasi¬ 
static  and  dynamic  experimental  results  at  different  strain  rates  at  room  temperature,  and 
high  &  low  temperature  experimental  results  at  one  strain  rate.  This  usually  provides  one 
with  a  set  of  initial  material  constants  which  can  be  used  as  input  into  a  computer 
software  which  utilizes  constrained  optimization  procedure  using  least  square  method  to 
obtain  a  set  of  more  refined  material  constants.  The  dynamic  data  before  being  input  in 
the  optimization  scheme  is  converted  into  an  equivalent  isothermal  data  by  calculating 
the  instantaneous  temperature  rise  in  the  material  at  any  strain  level  and  calculating  stress 
increment  for  that  temperature  rise  to  get  corresponding  isothermal  stress  response.  The 
material  constants  thus  determined  were  used  to  predict  a  strain  rate  jump  experiment, 
which  was  perfonned  separately  for  each  alloy.  This  experiment  was  not  included  in  the 
calculation  of  the  material  constants.  The  strain  rate  jump  experiment  included  at  least 
two  different  strain  rates  in  the  quasi-static  regime  followed  by  a  dynamic  regime  strain 
rate.  The  final  material  constants  for  the  three  alloys  are  given  below. 


A 

(MPa) 

B 

(MPa) 

Hi 

«e 

C 

m 

Ti-6A1-4V  Alloy  1 

1100 

857.5 

0.5455 

0.6086 

0.02204 

1.6236 

Ti-6A1-4V  Alloy  2 

988 

747.1 

0.5455 

0.3986 

0.02204 

1.2214 

Ti-6A1-4V  Alloy  3 

1069 

874.8 

0.5455 

0.4987 

0.02204 

1.3916 

Table  2.  KHL  model  material  constants  determined  for  the  three  Ti-6A1-4V  Alloys  used  in  the  investigation 

For  alloy  1,  five  different  strain-rate  experiments  were  perfonned  ranging  from 
10‘6  s'1  to  dynamic  strain  rate  of  3100  s'1.  The  material  behavior  was  found  to  be  strain 
rate  sensitive  and  has  a  non  linear  work  hardening  response  at  low  deformation  levels 
transitioning  to  an  almost  linear  response  at  higher  defonnation  levels.  Also,  there  is  a 
slight  reduction  in  work  hardening  rate  with  increasing  strain  rates  in  the  quasi-static 
loading  regime.  There  is  no  or  little  work  hardening  effect  seen  in  the  adiabatic  dynamic 
experiment  due  to  thermal  softening.  The  room  temperature  experiments  along  with  the 
model  correlations  can  be  seen  in  Fig  5.  The  observed  response  of  the  material  is  also  a 
non-linear  function  of  temperature.  The  increase  in  test  temperature  causes  a  drop  in  the 
flow  stress.  A  slight  decreasing  work  hardening  rate  behavior  can  be  noticed  as  the 
temperature  increases  and  also  with  an  increase  in  the  level  of  strain.  The  response  of  the 
alloy  at  various  temperatures,  along  with  the  model  correlations  for  different 
temperatures  can  be  seen  in  Fig  6. 


True  Stress  (MPa)  y>  True  Stress  (MPa) 


KHL  model  correlations  with  the  room  temperature  (296K)  response  of  alloy  1  at  different  strain 

rates. 


Fig.  6.  KHL  model  correlations  with  the  observed  response  of  alloy  1  at  different  temperatures  at  constant 

strain-rate  10"4  s"1 


Experiments  were  also  performed  at  high  temperatures  under  dynamic  loading 
conditions  on  alloy  3  and  the  response  was  modeled.  The  KHL  model  was  able  to 
successfully  predict  the  material  behavior  at  these  temperatures  and  strain  rates.  The 
point  to  note  is  that  these  set  of  experimental  observations  were  not  included  in 
calculating  the  material  constants.  The  experimental  observations  and  KHL  predictions 
are  shown  in  Fig  7.  Overall  KHL  model  is  able  to  simulate  different  loading  conditions  at 
different  strain  rates  and  temperatures  for  all  the  alloys  used  in  the  investigation. 


Fig.  7.  KHL  model  predictions  with  the  observed  response  of  alloy  3  at  high  strain  rates  and  different 

temperatures.  All  the  predictions  are  adiabatic. 


Executive  Summary  of  Achievements  III 

A  large  number  of  studies  have  been  carried  out  in  understanding  and  modeling  the 
response  of  these  titanium  alloys.  These  studies  were  focused  mainly  on  the  uniaxial 
response  of  the  alloys  at  different  strain  rates  and  temperatures,  and  assuming  isotropic 
behavior.  However,  it  is  well  known  that  the  stress  states  in  these  alloys  while  in  use  are 
definitely  more  complex  and  multiaxial  in  nature.  Hence  it  becomes  extremely  important 
to  understand  the  material  behavior  during  these  loading  conditions.  The  uniaxial 
material  responses  and  modeling  act  only  as  the  first  step  towards  understanding  the 


actual  more  complex  conditions  that  exist  and  how  they  affect  the  corresponding 
responses. 


Titanium  alloys  in  the  low  and  medium  temperature  regime  consist  mainly  of  the 
HCP  a  -phase  with  very  little  dispersed  (BCC)  (3  -phase  in  between  the  equiaxed  a 
grains.  These  hexagonal  closed  packed  (HCP)  metals,  especially  in  a  plate  geometry,  are 
known  to  display  plastic  anisotropy  and  a  strong  strength  differential  in  tension  and 
compression.  This  is  because  of  the  interaction  between  crystallographic  slip  and 
deformation  twinning  (Cazacu  et  ah,  2005).  Twinning  is  also  known  to  be  dependent  on 
the  sign  of  the  shear  stresses.  Also,  anisotropy  exists  in  these  alloys  mainly  due  to  the 
strong  crystallographic  texture  induced  due  to  the  rolling  process.  Therefore,  multiaxial 
modeling  of  these  HCP  titanium  alloys  cannot  be  performed  using  definitions  of  the 
effective  stress,  strain  and  strain  rate  based  on  the  isotropic  J?-flow  theory  as  these  alloys 
have  been  known  to  be  anisotropic  in  nature.  Similar  conclusions  were  found  by  other 
researchers  while  studying  the  multiaxial  response  of  a  -titanium. 

Anisotropy  has  been  studied  in  great  detail  by  many  researchers  over  the  past  many 
decades.  Hill  (1948)  was  the  first  to  quantitatively  determine  the  yield  surface  for 
anisotropic  materials.  Hill’s  anisotropic  model  was  primarily  based  on  the  generalization 
of  the  von-Mises  yield  criterion.  Although  this  model  has  been  useful  because  of  its 
simplicity  it  was  not  able  to  capture  the  observed  anisotropy  in  yield  stress  and  the  R- 
ratios  together  (Lademo  et  al.  1999).  Other  prominent  anisotropic  yield  criteria  are  by 
Hosford  (1979),  Barlat  and  co  workers  (Barlat  et  al.  1991,  1997,  2003),  Karafillis  and 
Boyce  (1993),  Yoon  et  al.  (2006),  etc.  These  criteria  dealt  with  the  modeling  of  cubic 
crystals  and  assume  same  yield  stress  in  tension  and  compression.  Some  anisotropic 
models  have  been  proposed  and  used  in  the  literature  for  HCP  poly-crystals  by  Tome  and 
Lebensohn  (2004),  and  Staroselsky  and  Anand  (2003).  Macroscopic  criteria  are  generally 
less  time  consuming,  such  as  the  orthotropic  criterion  by  Cazacu  et  al.  (2005)  which  is 
based  on  the  linear  transformation  of  the  Cauchy  stress  tensor  deviator. 

The  present  phase  of  investigation  includes  a  comprehensive  set  of  experiments  for 
understanding  the  response  of  this  alloy  and  consequently  modeling  the  observed 
responses  of  these  alloys  at  different  strain  rates.  It  is  an  extension  of  the  previous 
investigation  by  Khan  et  al.  (2004).  In  this  investigation,  the  high  strain  rate  torsion 
experiments  are  carried  out  using  the  modified  Kolsky  bar  apparatus  (Hartley  et  al., 
1985).  This  is  followed  by  compression  on  the  same  samples  again  at  high  strain  rates, 
albeit  at  different  operating  temperatures.  The  quasi-static  torsion  experiments  are  carried 
out  using  the  MTS  axial/  torsion  material  testing  system.  The  non-proportional  loading 
experiments  are  perfonned  using  the  channel  die  at  small  and  intennediate  strain-rate 
regime.  During  these  experiments  the  stress  state  is  designed  to  change  from  uniaxial  to 
biaxial  compression  and  the  corresponding  response  is  observed  and  recorded. 
Constitutive  modeling  of  these  observed  responses  is  performed  using  the  KHL  model. 
The  experimental  data  is  first  converted  into  effective  stress  and  effective  strain  based  on 
definitions  of  the  orthotropic  yield  criteria  proposed  by  Cazacu  et  al  (2005). 


In  order  to  describe  the  strength  differential  and  the  anisotropy  typically  shown  by 
hexagonally  closed  packed  (HCP)  metals  and  alloys,  Cazacu  et  al.  (2005),  proposed  a 
new  phenomenological  orthotropic  criterion  that  could  capture  the  response  of  the 
material  in  different  directions  accurately.  This  criterion  was  based  on  the  linear 
transformation  of  the  deviatoric  part  of  the  Cauchy  stress  tensor,  S  similar  to  previous 
studies  by  Barlat  and  co-workers  (Barlat  et  all 991;  Barlat  et  al.  1997;  etc)  and  had 
previously  been  used  by  many  researchers  (Lademo  et  al,  1999;  Lademo  et  al,  2002; 
Abedrabbo  et  al.,  2006,  etc)  successfully.  The  transformed  tensor  Z  is  defined  as: 

Z  =  C[s]  (3) 

where  C  is  a  4  th  order  tensor.  The  ortho  tropic  yield  criterion  is  given  as: 

(II,  |  -  *.I, )“  +  QX2|  -  kX2 )“  +  (p,  |  -  i.X, )"  =  F  (4) 

The  exponent  a,  in  the  yield  criterion  is  considered  to  be  a  positive  integer  (in  our 
case  a= 2).  The  parameter  k  is  a  material  constant,  which  for  a  fixed  a,  is  detennined  by 
the  ratio  of  the  strength  in  tension  and  compression.  The  function  F,  gives  the  size  of  the 
yield  locus.  F  can  be  defined  either  as  a  constant  or  a  function  of  the  total  plastic  strain 
(isotropic  hardening). 

Experimental  results  from  uniaxial  tension  and  compression  yield  stresses,  in 
rolling  and  transverse  to  rolling  directions,  uniaxial  compression  in  the  45°  to  the  rolling 
direction  and  R-ratios  in  the  rolling,  and  transverse  to  rolling  directions  for  the 
compression  case  are  used.  These  material  constants  are  calculated  using  least  square 
optimization  technique  in  MATLAB.  The  value  of  the  material  constant  k,  used  in  this 
investigation  is  calculated  from  the  ratio  of  the  yield  in  tension  and  compression  in  the 
transverse  to  rolling  direction.  The  experimental  yield  stresses  for  rolling  and  transverse 
to  rolling  as  well  as  the  yield  surface  generated  from  the  orthotropic  criteria  are  shown  in 
Fig. 8.  The  corresponding  experimental  R-ratios  and  the  yield  stresses  at  45°  to  the  rolling 
direction  and  in  shear  are  shown  along  with  the  correlations  from  the  model  in  Table  3. 


Fig.  8.  Yield  surface  generated  using  the  orthotropic  yield  criteria  from  Cazacu  et  al.  (2005)  for  Ti-6A1-4V 

titanium  alloy 


Yield  Stress 

aT0 

(TT90 

aCo 

„c 
o  90 

O  45 

Oxv 

rc 
r  0 

rc 

r  90 

Experimental  Values 

138.61 

130.9 

157.67 

135.44 

144.01 

74.76 

0.685 

0.507 

Normalized  Exp.  Values 

1.02 

0.97 

1.16 

1 

1.06 

0.55 

0.685 

0.507 

Predicted  Values 

1.04 

0.94 

1.14 

1.03 

1.04 

0.58 

0.69 

0.501 

Table  3.  Experimental  and  predicted  values  of  the  yield  stresses  and  R-ratios  under  different  loading 
conditions 


It  can  be  seen  from  the  Fig.  8  and  Table  3  that  the  orthotropic  model  works 
reasonably  well  for  HCP  metals  in  the  development  of  the  yield  surface  and  that  the 
values  of  R-ratios  (which  incorporates  the  flow  rule)  correlated  by  this  model  are 
extremely  close  to  the  experimental  data.  This  implies  that  it  is  reasonable  to  use 
effective  stress,  strain  and  strain  rate  definitions  using  this  ortho  tropic  model. 


The  experimental  results  from  quasi-static  mo  no  tonic  torsion,  quasi-static  jump 
torsion  and  dynamic  torsion  experiments  were  converted  into  effective  stress  ( <7e )  using 
the  definition  of  the  yield  function  (Eqn.  4),  so  that  the  complexity  of  the  stress  state 
reduces  to  a  single  scalar  quantity.  In  a  similar  manner  the  effective  strain,  ( £e ),  and 

strain  rates,  ( £e ),  are  also  calculated  from  the  experimental  results  by  considering  the 

principle  of  incremental  work  per  unit  volume  (Khan  and  Huang,  1995).  The 
experimental  results  were  then  plotted  against  the  model  predictions  based  on  these 
definitions. 


dWp  =  (Jjdsjj  =  oid£r’  =  (7  ( cl  £p  (5) 

Fig.  9  shows  the  experimentally  observed  multiaxial  response  of  the  alloy  along 
with  the  model  predictions  (the  material  constants  were  not  changed  from  the  published 
values,  Khan  et  ah,  2004).  The  KHL  model  is  able  to  predict  the  response  of  these  quasi 
static  as  well  as  the  dynamic  torsion  experiments  reasonably  well.  The  KHL  predictions 
are  made  assuming  adiabatic  deformations  for  the  1.732  s'1  (shear  strain  rate)  as  well  as 
the  dynamic  experimental  results  as  an  increase  in  the  temperature  were  noticed  during 
both  experiments. 


Fig.  9.  Model  predictions  with  the  observed  shear  response  of  the  titanium  alloy  converted  to  effective 
stress  and  strain  during  monotonic,  strain-rate  jump  and  dynamic  torsion  experiments  at  room  temperature 
(296K).  Dynamic  and  0.989  s"1  effective  strain  rate  predictions  are  based  on  assuming  adiabatic 

deformation. 


Fig.  10  shows  the  measured  response  of  the  alloy  under  dynamic  torsion  at  room 
temperature,  followed  by  dynamic  compression  at  higher  temperature  (588  K).  The  KHL 
model  predictions  were  converted  for  adiabatic  case  as  was  done  in  the  previous  case. 
The  predictions  are  again  in  very  good  agreement  with  the  actual  observed  material 
response  for  both  types  of  dynamic  loading  conditions  showing  the  capability  of  the 
model  in  capturing  the  observed  behavior  of  the  alloy  under  complex  loading  conditions 
when  used  in  conjunction  with  the  anisotropic  criterion. 


Fig.  10.  Dynamic  shear  response  at  room  temperature,  followed  by  dynamic  compression  response  at 
588°K  of  the  titanium  alloy,  and  predictions  using  KHL  model  with  the  material  constants  determined  from 
uniaxial  compression  experimental  data. 


The  true  stress  and  strain  response  of  the  alloy  under  non-proportional  loading 
conditions,  uni-axial  to  bi-axial  compression,  at  10'1  s’1  strain  rate  in  the  loading  direction 
is  shown  in  Fig.  11.  The  rise  in  the  stress  in  the  loading  direction  is  clearly  visible  as  the 
specimen  becomes  constrained  once  a  certain  predetermined  level  of  deformation  in  the 
loading  direction  is  achieved.  Consequently,  there  is  a  corresponding  rise  in  the  stress  in 
the  constrained  direction  until  the  end  of  the  experiment.  Fig.  12  shows  the  effective  stress 
and  effective  strain  responses  for  different  strain  rates  in  the  loading  directions  along  with 
the  model  predictions.  Effective  stresses,  strains  and  strain  rates  are  calculated  based  on 
the  orthotropic  yield  criterion  of  Cazacu  et  al  (2005)  and  the  principle  of  incremental 
work  per  unit  volume.  The  stresses  and  strains  are  then  plotted  with  the  model 
predictions.  Again,  the  published  values  of  constants  [Khan  et  al  (2004)]  are  used  and 
are  not  modified  using  the  present  experimental  results.  The  model  prediction  for  the  case 


of  10'1  s'1  strain  rate  in  the  loading  direction  is  converted  assuming  adiabatic  deformation, 
as  there  was  a  significant  increase  in  the  temperature  in  the  specimen  noticed  during  the 
deformation.  The  model  predictions  are  in  excellent  agreement  with  the  experimentally 
observed  responses  for  all  strain  rates.  The  strain  rates  used  in  the  model  calculations  in 
all  non  proportional  loading  experiments  are  effective  strain  rates,  and  are  computed 
using  the  three  strain  histories  along  with  the  principle  of  incremental  work  per  unit 
volume. 


Fig.  11.  True  stress-strains  for  the  non-proportional  loading  channel  die  experiment  at  10"'  s"1  strain  rate  in 
the  loading  direction.  Note  the  stress  starts  to  increase  in  the  direction  2  (constrained  direction)  from  zero 
once  the  specimen  comes  in  contact  with  the  die  walls  at  approximately  12%  strain  in  the  loading  direction. 


Fig.  12.  Quasi  static  response  of  the  titanium  alloy  under  non  proportional  biaxial  loading  at  room 
temperature  converted  to  effective  stress  and  strain  based  on  an  orthotropic  yield  criteria  and  predictions 
using  KHL  model  with  the  material  constants  determined  from  uniaxial  compression  experimental  data. 

Using  the  set  of  model  constants  for  this  alloy  which  were  calculated  from  the 
uniaxial  responses  at  different  strain  rates  and  temperatures  in  an  earlier  paper  (Khan  et 
al.  2004),  the  capability  of  the  model  to  predict  accurately  the  observed  new  responses 
under  the  complex  multiaxial  loading  conditions  was  demonstrated.  The  model  was  able 
to  predict  the  material  response  under  pure  torsion,  and  torsion  followed  by  compression 
extremely  well  using  the  effective  stress,  strain  and  strain  rate  definitions  from  a 
orthotropic  yield  criterion  and  principle  of  incremental  work  per  unit  volume.  In  the  non 
proportional  channel  die  experiments,  the  model  was  again  able  to  predict  the  observed 
material  responses  under  non-proportional  loading  conditions  reasonably  well. 
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Abstract 

The  results  from  a  systematic  study  of  the  response  of  a  Ti-6A1^1V  alloy  under  quasi¬ 
static  and  dynamic  loading,  at  different  strain  rates  and  temperatures,  are  presented.  The 
correlations  and  predictions  using  modified  Khan-Huang-Liang  (KHL)  viscoplastic 
constitutive  model  are  compared  with  those  from  Johnson-Cook  (JC)  model  and 
experimental  observations  for  this  strain  rate  and  temperature-dependent  material. 
Overall,  KHL  model  correlations  and  predictions  are  shown  to  be  much  closer  to  the 
observed  responses,  than  the  corresponding  JC  model  predictions  and  correlations. 
Similar  trend  has  been  demonstrated  for  other  titanium  alloys  using  published 
experimental  data  [Mech.  Mater.  33(8)  (2001)  425;  J.  Mech.  Phys.  Solids  47(5)  (1999) 
1157]. 
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1.  Introduction 

Dynamic  deformation  has  been  of  interest  not  only  in  impact  and  penetration  related 
problems  (  Gilat  &  Cheng,  2002;  Khan  et  al.  2002;  Khan  &  Lopez-Pamies,  2002;  Bjerke 
et  al.  2002)  but  also  in  high  speed  machining  (Molinari  et  al.  2002;  Burns  &  Davies  2002; 
Titanium  alloys  have  been  studied  by  several  investigators  because  of  their  use  in  aero¬ 
engine,  gas  turbines  and  other  applications  due  to  their  high  strength  to  weight  ratio, 
ductility,  and  ability  to  withstand  high  temperatures  and  resist  corrosion.  Examples  are 
Macdougall  and  Harding  (1999),  Meyers  et  al.  (1994),  Chichili  et  al.  (1998),  Cheng  and 
Nemat-Nasser  (2000).  The  first  investigation  mentioned  above  was  on  Ti-6A1-4V  alloy, 


while  others  were  on  commercially  available  pure  alpha-titanium  except  in  case  of  Cheng 
and  Nemat-Nasser  (2000),  in  which  the  microstructure,  or  phase,  was  not  specified.  The 
development  of  relatively  economical  Ti-6A1-4V  alloy,  with  resulting  high  oxygen 
content,  has  sparked  interest  in  its  possible  use  in  lightweight  tanks  (Montgomery  and 
Wells,  2001);  the  conventional,  more  expensive  Ti-6A1-4V  alloy  has  been  used  primarily 
in  aerospace  components.  The  potential  applications  in  armor,  including  ceramic  tiles 
encapsulated  in  titanium  alloys,  have  motivated  several  studies  (Gray,  1997;  Follansbee 
and  Gray,  1989;  Lesuer,  2000;  Nemat-Nasser  et  ah,  2001;  Majorell  et  ah,  2002). 

In  case  of  a-titanium,  it  has  been  demonstrated  that  the  deformation  mechanisms 
include  glide  systems  with  a-type  dislocations  in  the  h.c.p.  structure,  as  well  as,  twinning 
shear  which  contributed  to  the  overall  strain  (Meyers  et  ah,  1994;  Song  and  Gray,  1995; 
Chichili  et  ah,  1998).  In  the  only  quantitative  study  of  twinning,  Chichili  et  ah  (1998) 
have  shown  that  twin  number  density  increases  drastically  with  increase  in  strain  rate. 
They  used  length-to-diameter  ratio  of  1.6  for  quasi-static  experiments.  The  rate  of 
increase  of  this  density  decreased  with  increase  in  strain  at  dynamic  strain-rates  of  10  s'  , 
while  this  rate  increased  with  deformation  in  the  case  of  quasistatic  loading  (10'5  s'1). 
Their  results  further  demonstrated  that  loading  at  a  particular  temperature  (and/or  strain 
rate),  unloading,  and  reloading  at  another  temperature  (and/or  strain  rate),  resulted  in  a 
stress-strain  curve  which  was  significantly  different  than  if  the  specimen  was  loaded  at 
the  latter  temperature  (and/or  strain  rate)  right  from  the  beginning.  This  study  was 
performed  over  a  wide  range  of  strain  rates  (10'  -10  s'  )  but  over  a  very  limited  range  of 
temperatures  (77-298  K).  Further,  they  did  not  provide  the  geometry  of  the  specimens 
used  in  their  split  Flopkinson  bar  measurements. 

The  Ti-6A1-4V  alloy  consists  of  hep  a-grains,  with  a  dispersion  of  stabilized  bee  P 
phase  around  grain  boundaries  at  room  temperature,  a  phase  transfonns  to  P  phase 
starting  at  873  K  (1110  F);  above  1268  K  (1825  F),  the  entire  microstructure  consists  of 
equiaxed  P  grains  (Majorell  et  ah,  2002).  This  alloy,  in  addition  to  H,  V,  and  Ti,  contains 
O,  Fe,  Mo,  C,  Si,  and  Mn.  Oxygen,  nitrogen,  and  carbon  contents  are  a  stabilizers 
(Conrad  et  ah,  1975),  while  vanadium,  iron  and  molybdenum  are  P  stabilizers.  Conrad  et 
ah  (1975)  proposed  an  equivalent  oxygen  content  (Oeq  =  O+2N+0.75C);  this  equivalent 

impurity  concentration  gives  the  effect  of  dislocations-impurity  interaction  on  the  yield 
strength  of  the  material.  Investigation  by  Majorell  et  ah  (2002)  was  on  an  untextured  and 
a  textured  Ti-6A1-4V  alloy  rod  that  was  manufactured  by  Allvac.  This  study  was  over  a 
strain-rate  range  of  10'  to  10  s'  and  a  temperature  range  of  650-1345  K  (710-1970  F). 
They  did  not  observe  any  dynamic  strain-aging  at  any  of  the  temperatures  or  strain  rates 
investigated.  Further,  athennal  stress,  i.e.,  temperature  insensitive  response,  was  found  at 
approximately  1255  K  (1800  F);  i.e.  at  a  temperature  when  the  material  contained  almost 
100  %  P  phase.  They  did  not  specify  specimen  geometries  used  in  their  investigation.  The 
equivalent  oxygen  content  was  0.206  %. 

In  a  study  by  Follansbee  and  Gray  (1989),  on  a  Ti-6A1-4V  alloy,  with  an  equivalent 
oxygen  content  of  0.23%  (actual  oxygen  was  0.18%),  the  measurements  were  restricted 
to  three  temperatures  with  a  range  between  76  and  295  K  and  at  only  two  strain  rates 


3  1 

(10"  and  approx.  3000  s'  ).  Their  specimens  on  as-received  and  two  heat-treated  versions 
were  made  from  a  13.8  mm  thick  plate;  specimens  were  4.8  mm  in  diameter  and  5.2  mm 
long  in  case  of  quasistatic  loading  experiments,  and  6.4  mm  diameter  with  a  length  of  5.1 
mm  in  a  12.7  mm  diameter  split  Hopkinson  bar  experiment.  This  investigation  did  not 
include  responses  over  a  wide  range  of  temperatures.  Nemat-Nasser  et  al.  (2001)  study 
was  on  a  commercial  and  two  hot  isostatically  pressed  Ti-6A1-4V  alloys.  The  equivalent 
oxygen  in  the  commercially  pure  alloy  was  0.22  %.  The  geometry  of  the  specimens  was 
not  specified  and  was  presumably  5  mm  in  diameter  and  5  mm  long,  for  quasistatic  and 
dynamic  experiments.  The  diameters  of  the  split  Hopkinson  bars  (perhaps  12.7  mm)  were 
not  specified  either  in  the  paper.  The  measurements  and  assumptions  in  this  investigation 
have  raised  several  questions.  First,  they  perform  experiments  on  specimens  from  77  K  to 
998  K  in  the  dynamic  strain  rate  regime  of  2000  to  4000  s'1.  As  discussed  earlier,  a-phase 
starts  transforming  to  [3-phase  at  873  K.  Thus,  they  performed  experiments  on 
presumably  different  materials  in  their  range  of  temperatures,  i.e.,  with  different  amounts 
of  a  and  [3  phases.  Second,  it  is  clear  from  the  figures  in  their  paper  that  at  least  some 
experiments  had  barreling  due  to  insufficient  lubrication.  Third,  they  combine  initial 
portions  of  multiple  stages  loading-unloading-reloading  experiments  using  split 
Hopkinson  bar  technique,  to  obtain  “isothermal”  response,  even  when  it  is  well  known 
that  initial  portion(s)  of  the  stress-strain  curve  is  not  accurate  as  force  equilibrium  in  the 
specimen  is  not  reached  until  elastic  and  plastic  waves  are  reflected  several  times  back 
and  forth  in  the  specimen.  Using  these  loading-unloading  and  reloading  experiments  at 
different  temperatures,  they  suggested  that  all  (i.e.  100  %)  of  plastic  work  done  was 
converted  to  heat.  This  suggestion  is  in  direct  contradiction  to  the  measurements  of 
Mason  et  al.  (1994),  Liao  and  Duffy  (1998),  Macdougall  and  Harding  (1999)  and  Rosakis 
et  al.  (2000).  As  mentioned  earlier,  Chichili  et  al.  (1998)  and  Follansbee  and  Gray  (1989) 
have  shown  that  the  response  of  a  material  is  different  upon  reloading  to  a  different 
temperature  (same  strain  rate),  presumably  due  to  history  effects. 

Similar  inadequacies  and  contradictions  are  found  in  the  published  modeling  efforts. 
The  constitutive  models  used  for  high  strain  rate  applications  can  be  classified  in  two 
categories;  the  purely  phenomenological  ones,  e.g.  Johnson-Cook  (J-C)  (Johnson  and 
Cook,  1983)  and  Khan-Huang-Liang  (KHL)  models  (Khan  and  Huang,  1992;  Khan  and 
Liang,  1999;  Khan  and  Zhang,  2000,  2001)  and  so  called,  “physically  based  models”  e.g. 
the  ones  by  Zerilli  and  Armstrong  (1987),  Mecking  and  Kocks  (1981),  etc.,  that  were 
used  frequently  by  Follansbee  and  Gray  (1989),  and  in  a  modified  form  by  Nemat-Nasser 
et  al.  (2001).  The  latter  group  discusses  the  mechanisms  of  plastic  defonnation,  mainly 
dislocations.  However,  the  material  constants  are  determined  by  not  measuring  any 
deformation  mechanism  related  quantity,  but  by  choosing  constants  to  “fit”  the  uniaxial 
stress-strain  curves  at  different  strain  rates  and  temperatures,  just  like  the  purely 
phenomenological  models.  In  order  to  “fit”  their  model  to  the  uniaxial  stress-strain  curves 
of  two  fairly  close  Ti-6A1-4V  alloys,  Nemat-Nasser  et  al.  (2001)  used  athennal  stress  of 
685  MPa  while  Follansbee  and  Gray  (1989)  used  100  MPa;  Follansbee  and  Gray  (1989) 
used  a  value  of  10  MPa  to  analyze  results  of  Paton  et  al.  It  is  disappointing  to  note  that 
almost  six  hundred  percent  difference  in  this  value  by  these  two  set  of  researchers  using 
models  that  are  presumably  based  on  the  same  “physics”  for  very  similar  materials.  In  the 
original  Mecking  and  Kocks  model,  as  used  by  Follansbee  and  Gray  (1989),  a  highly 


undesirable  extrapolation  to  0  K  is  necessary.  In  addition,  the  model  has  numerous 
material  constants;  it  requires  huge  amount  of  resources  and  time  to  determine  these 
constants.  In  the  simplification  or  modification  of  this  model  by  Nemat-Nasser  et  al. 
(2001),  experiments  are  needed  at  very  high  temperatures  where  proper  lubrication  is 
extremely  difficult,  if  not  impossible,  to  attain,  as  is  clearly  evident  in  their  paper.  Since 
both  categories  of  models,  phenomenological  and  “so  called  physically  based”,  detennine 
material  constants  by  “fitting”  to  the  stress-strain  responses  at  different  temperatures  and 
strain  rates,  the  advantage  of  one  over  the  other  is  merely  number  of  material  constants  in 
each  model.  Any  model  with  semi-infinite  number  of  material  constants  will  be  able  to 
approximate  observed  responses.  Mecking  and  Kocks  model,  as  used  by  Follansbee  and 
Gray  (1989)  has  23  constants,  while  these  constants  range  from  12  to  8  in  case  of  Cheng 
and  Nemat-Nasser  (2000)  and  Nemat-Nasser  et  al.  (2001),  respectively,  depending  on 
whether  they  include  modeling  of  dynamic  strain  aging  or  not.  Johnson-Cook  and  Khan- 
Huang-Liang  models  have  5  and  6  constants,  respectively.  Johnson-Cook  and  Khan- 
Huang-Liang  models  are  used  in  this  investigation  due  to  their  advantage  of  fewer 
constants  and  their  ability  to  model  the  observed  material  response  as  closely  as  with 
models  with  many  more  constants. 

Modeling  Approach 


The  material  constants  for  both,  the  modified  KHL  (Khan-Huang-Liang)  and  JC 
(Johnson-Cook)  models  were  determined.  Using  these  material  constants,  correlations 
were  obtained  and  compared  to  experimental  results.  The  modified  Khan-Huang-Liang 
(KHL)  model  is  as  follows. 


<j  =  [A  +  B{  1- 


In  P  P  T* 

In  D%’  V  Tm 


(1) 


where  o  is  the  stress  and  £p  is  the  plastic  strain.  Tm,  T ,  Tr  are  melting,  current,  and 

reference  temperatures,  respectively.  Dp=  106  s'1  (arbitrarily  chosen  upper  bound  strain 

rate)  and  £=\  s'1  (reference  strain  rate  at  which  some  material  constants  are  determined). 
£  is  the  current  strain  rate.  A,B,nvn0,C  and  m  are  material  constants.  For  Ti-6A1-4V 

alloys,  the  melting  temperature  is  1933  K  (Nemat-Nasser  et  al.,  2001).  A  major  feature  of 
this  model,  unlike  the  Johnson-Cook  (JC)  model,  is  that  decreasing  work-hardening  with 
increasing  strain  rate  can  be  accommodated  through  the  material  constant  n] .  Previously, 
this  model  had  a  temperature  dependence  term  (similar  to  JC  model)  as: 


l-( 


T-T 
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However,  this  tenn  can  not  accommodate  the  case  when  the  current  temperature  is  lower 
than  the  reference  temperature  as  then  the  above  term  within  parenthesis  becomes  a 


negative  number  raised  to  the  power  m  .  Therefore  a  modification  was  made  as  given  in 
Equation  (1).  The  Johnson-Cook  (JC)  model  is  given  by  the  following  equation. 


(j  =  [A  +  B(£p)n°\l  +  C\n—  ) 
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In  this  case,  various  terms  have  the  same  definitions  as  Equation  (1). 

To  determine  the  material  constants  for  the  model,  uniaxial  loading  test  results  at 
different  strain  rates  and  temperatures  were  first  used  to  obtain  an  initial  set  of  values  for 
the  material  constants,  as  described  in  Khan  and  Liang  (1999).  In  order  to  get  a  refined 
set  of  material  constants,  the  least  squares  and  constrained  optimization  procedures  were 
used  to  minimize  the  difference  between  correlated  and  actual  data  using  MATLAB.  The 
experimental  results  at  low  strain  rates  were  first  used  to  get  a  set  of  material  constants 
using  the  least  squares  method.  Then,  higher  strain  rate  data,  measured  under  adiabatic 
condition,  were  correlated  with  the  assumption  that  90  %  plastic  work  was  dissipated  to 
heat,  to  get  a  better  set  of  material  constants  for  both  low  and  high  strain  rate  material 
response  using  a  constrained  optimization  procedure.  This  was  iteratively  performed  until 
the  constants  converged.  As  for  the  dynamic  data,  thennal  softening  from  the  adiabatic 
deformation  was  effectively  considered  by  converting  the  increment  of  temperature  from 
the  stress-strain  curve  using  following  equation: 


B  r£i 

AT  =  — f '  cr{ep)dep  (4) 

PCP  Jo 

where  f3,  p,  Cp  are  the  fraction  of  heat  dissipation  caused  by  the  plastic  deformation, 

mass  density  and  specific  heat  at  constant  pressure,  respectively.  In  fact,  there  are  many 
arguments  on  the  conversion  of  plastic  work  to  heat  during  high  strain  rate  defonnation. 
Taylor  and  Quinny  (1934)  measured  the  fraction  of  heat  conversion  and  found  that  it  is  of 
the  order  of  0.9  with  copper  specimens.  Mason  et  al.  (1993)  measured  it  with  dynamic 
experiments  and  demonstrated  that  it  depended  substantially  on  strain  and  strain  rate. 
Kapoor  and  Nemat-Nasser  (1998)  measured  the  energy  converted  to  heat  using  an  infra¬ 
red  method  for  Ta-2.5  %  W  alloy,  commercially  pure  Ti,  1018  steel,  6061  Al  and  OFHC 
Cu.  They  argued  that  the  infra-red  measurement  generally  underestimated  the  conversion 
ratio  (70%  of  conversion  of  work  for  Ta-2.5  %  W,  for  example),  they  cooled  down  the 
specimen  to  the  initial  room  temperature  at  certain  strain  and  heated  it  to  the  temperature 
based  on  100%  heat  conversion  and  then  deformed  the  specimen  at  the  same  strain  rate  to 
check  if  the  adiabatic  curve  was  obtained.  On  this  basis,  they  suggested  that  nearly  all  of 
the  plastic  work  was  converted  into  heat  (i.e.,  /)  =  1 )  within  experimental  error, 
concluding  that  the  infra-red  detection  system  recorded  a  lower  (surface)  temperature 
than  the  actual  temperature  of  the  specimen.  Later,  Nemat-Nasser  et  al.  (2001)  argued  for 
this  total  heat  conversion  by  measuring  the  thermomechanical  response  of  a  Ti-6A1-4V 
alloy  (commercial  Ti64).  Macdougall  and  Harding  (1999)  calculated  the  proportion  of 
work  converted  to  heat  based  on  their  torsional  dynamic  test  with  an  infrared  technique 


for  Ti-6A1-4V  specimens  and  concluded  that  (3  varied  with  increasing  plastic  strain  from 
approximately  0.2  to  approaching  0.7.  Rittel  (1999)  and  Rosakis  et  al.  (2000)  also 
indicated  that  (3  was  a  history-dependent  quantity  and  depended  strongly  on  both  strain 
and  strain  rate  for  various  engineering  materials.  There  seems  to  be  no  agreement  on  a 
single  value  of /3  that  works  for  variety  of  engineering  materials;  however  most 
researchers  still  use  0.9,  e.g.  Lennon  &  Ramesh  (2004).  In  this  study,  (3  =  0.9  was 
assumed. 

For  Ti-6A1-4V,  the  density  is  p  =  4428  kg/m3  (Lesuer,  2000).  The  heat  capacity  was 
expressed  as  a  function  of  the  temperature,  (Military  Handbook,  1998); 


C=  559.77  -  0.1473T  +  0.00042949T2  J/(kg-K)  (278  K<  T  <  1144  K)  (5) 

Finally,  using  these  material  constants  for  both  models,  predictions  were  made  and 
compared  to  a  strain-rate-jump  experiment  (10“5  to  10'1  to  1700  s'1). 

Material  constants  for  both  models  were  detennined  for  two  additional  Ti-6A1-4V 
alloys  (Nemat-Nasser  et  al.,  2001;  Macdougall  and  Harding,  1999).  The  first  one  was 
commercially  available  from  Protech  Metals  (Nemat-Nasser  et  al.,  2001),  and  the  second 
one  was  used  in  torsion  and  tension  tests  by  Macdougall  and  Harding,  (1999),  where  the 
exact  composition  of  the  Ti-6A1-4V  was  not  given. 

Experimental  Details 


3.1  Material 

The  chemical  composition  of  Ti-6A1-4V  is  shown  in  Table  1.  This  Ti-6A1-4V  alloy  is 
ductile  and  has  a  hexagonal  close  packed  crystalline  structure.  The  composition  is 
slightly  differed  from  the  commercial  Ti  alloy  that  was  used  by  Nemat-Nasser  et  al 
(2001),  and  is  presented  in  Table  2. 


3.2  Experimental  Procedure 


3.2.1  Room  Temperature  Compression  Experiments 

Quasistatic  compression  experiments  were  conducted  using  an  MTS  servo  hydraulic 
axial/torsional  material  testing  machine  at  strain  rates  of  10'5  s'1  to  1  s'1,  using  cylindrical 
specimens.  The  specimens  were  machined  from  a  plate  (27.9  mm  thick)  of  Ti-6A1-4V  to 

19.1  mm  in  length  and  12.7  mm  in  diameter.  The  length  dimension  was  in  the  thickness 
direction  of  the  plate.  It  should  be  noted  that  the  experiment  at  a  strain  rate  of  1  s'1  was 
conducted  in  increments  of  5  %  strain,  with  30  minutes  between  loading  and  reloading. 
This  method  allowed  the  specimen  sufficient  time  to  cool,  ensuring  an  isothermal 


response.  For  room  temperature  compression  experiments,  KFEL-2-120-C1  high 
elongation  uniaxial  strain  gages,  manufactured  by  Kyowa  Ltd.  (Japan),  were  mounted 
diametrically  opposite  to  each  other  at  the  center  of  each  specimen  using  Micro- 
Measurements  Inc.  adhesive  AE-15.  The  MTS  transducers  directly  supplied  load  and 
LVDT  output  data  to  the  TestStar  software.  The  interface  between  the  load  platens  and 
the  test  specimens  for  each  of  the  compression  experiments  was  lubricated  with  Teflon 
sheets  and  Dow  Corning  high  vacuum  grease  in  order  to  achieve  homogeneous 
deformation  and  thus  avoid  barreling  of  the  specimen,  maintaining  a  uniform,  uniaxial 
stress  state. 


3.2.2  Compression  Experiments  at  Different  Temperatures 

Quasistatic  experiments  at  different  temperature  were  conducted  at  233,  296,  422, 
588,  and  755  K  with  a  strain  rate  of  10'  s'  .  The  LVDT  displacement  supplied  stroke  data 
that  included  the  total  defonnation  of  the  MTS  actuators,  the  loading  platens,  the  grease 
lubricant,  and  the  specimen.  The  interface  between  the  load  platens  and  the  test 
specimens  for  each  of  the  compression  experiments  was  lubricated  with  high  or  low 
temperature  grease,  depending  on  the  experiment  temperature,  in  order  to  achieve 
uniformity  in  defonnation  and  avoid  barreling  of  the  specimen,  maintaining  a  unifonn, 
uniaxial  stress  state.  The  stroke  was  corrected  by  running  a  blank  experiment  (without 
specimen). 


3.2.3  Dynamic  Compression  Experiments 

The  Split  Hopkinson  Bar  technique  was  used  to  attain  dynamic  measurements.  The 
interfaces  between  the  incident  and  transmitted  bars  and  the  test  specimen  were  carefully 
lubricated  with  grease,  in  order  to  increase  uniformity  in  deformation.  The  initial 
thickness  and  diameter  of  the  specimen  were  5.1  (in  the  direction  of  thickness  of  the 
plate)  and  10.2  mm,  respectively.  The  pressure  and  striker  bars  used  for  the  experiment 
were  made  of  hardened  VascoMax  C350  steel  and  the  diameters  were  12.7  mm. 


3.2.4  Room  Temperature  Strain-Rate-Jump-Compression  Experiments 

In  order  to  conduct  the  strain-rate -jump  experiment,  the  specimen  was  loaded  at  an 
initial  strain  rate  of  10'5  s'1.  The  strain  rate  was  then  increased  to  10"1  s'1  around  6.3  % 
strain.  Upon  unloading  the  specimen  after  quasistatic  compression,  around  14%  strain,  it 
was  machined  into  two  dynamic  compression  specimens  with  diameter  10.2  mm  and 
thickness  5.1  mm.  The  split  Hopkinson  pressure  bar  technique  was  then  used  to  obtain 
the  response  of  the  material  in  the  third  stage  of  loading  at  a  strain  rate  of  1700  s'1. 


Results  and  Discussion 


4.1  Experimental  Results 

Fig.  1  includes  the  measured  true  stress-true  plastic  strain  curves  shown  with  symbols 
(the  plastic  strain  was  converted  from  the  total  strain  measurement)  at  different  strain 
rates  and  at  room  temperature.  As  seen  in  Fig.  1,  the  material  response  at  various  strain 
rates  indicates  that  the  material  is  strain  rate  dependent.  Fig.  2  gives  the  measured  true 
stress  -  true  plastic  strain  curves  (also  shown  with  symbols)  at  different  temperatures  and 
at  a  strain  rate  of  10"  s'  .  Figures  1  to  3  include  results  using  the  Alloy  1  (see  Table  1). 
Figure  4  contains  results  of  Nemat-Nasser  et  al.  for  Alloy  2,  while  Figures  5  &  6 
represents  results  of  Macdougall  &  Harding  (1999). 


4.2  Determination  of  Material  Constants 

In  order  to  determine  the  material  constants,  all  data  were  included  in  the  least 
squares  method  and  constrained  optimization  procedure.  Isothermal  dynamic  response 
was  obtained,  and  included  in  the  optimization  procedure,  from  adiabatic  measurements 
using  /?  =  0.9.  In  case  of  JC  model,  the  material  constants  were  detennined  by  two 
methods,  the  published  non-optimized  approach  and  also  the  optimized  method  used  for 
KHL  method.  Tables  3  and  4  include  material  constants  for  KHL  and  JC  models, 
determined  for  Ti-6A1-4V  alloy  and  for  two  alloys  from  the  studies  of  the  data  extracted 
from  Nemat-Nasser  et  al.  (2001)  and  Macdougall  and  Harding  (1999). 

4.2.1  Correlations  and  Predictions  with  Ti-6Al-4V  (Alloy  1) 

The  quasistatic  and  dynamic  measurements  at  a  temperature  of  296  K  and  at  various 
strain  rates,  along  with  correlations  by  KHL  and  JC  models  are  shown  in  Fig.  1.  Although 
the  correlations  with  both  models  are  reasonably  good  for  the  quasistatic  case,  the  KHL 
model  correlates  better  with  the  dynamic  response,  while  the  JC  model  differs 
substantially  from  the  measured  response  in  terms  of  initial  yield  stress  and  work 
hardening  rate  at  this  high  strain  rate  defonnation.  It  was  observed  that  even  with  the 
optimized  material  constants,  JC  model  did  not  correlate  well  with  the  dynamic  response. 
This  is  due  to  the  fact  that  JC  model  fails  to  correlate  the  decreasing  work  hardening  rate 
as  the  strain  rate  increases.  This  is  described  elsewhere  in  detail  (Liang  and  Khan,  1999), 
in  which  it  was  shown  that  in  the  JC  model,  the  work  hardening  rate  at  certain  strain 
would  increase  when  strain  rate  increases.  That  is,  the  JC  model  is  not  appropriate  for 
modeling  any  material  where  the  work  hardening  rate  decreases  with  increasing  strain 
rate,  such  as  tantalum,  (Chen  and  Gray,  1995).  With  Ti-6A1-4V,  too,  it  can  be  observed 
from  the  material  behavior  at  different  strain  rates  that  the  work  hardening  rate  decreases 
as  the  strain  rates  increases. 

In  Fig.  2,  quasistatic  measurements  (at  a  strain  rate  of  10"  s'  )  at  various 
temperatures,  with  correlations  using  KHL  and  JC  models  are  shown,  respectively. 
Material  constants  for  JC  model  determined  using  published  conventional  methods,  i.e., 


without  the  optimization,  led  to  slight  deviation  from  the  measured  response.  Overall,  the 
KHL  model  gave  somewhat  better  correlation  than  the  JC  model,  except  for  422  K.  The 
JC  model  has  a  major  drawback  that  it  cannot  be  used  at  temperatures  lower  than  the 
reference  temperature  at  which  material  constants  are  determined.  Note  the  absence  of 
correlation  at  233  K  in  JC  model  as,  (T  -  Tref)/(Tm  -  Tref)  becomes  negative. 

A  strain-rate-jump  experiment  was  performed  at  296  K  with  strain  rate  from  10'5  s'1 
to  10"1  s'1,  followed  by  an  experiment  at  a  strain  rate  of  1700  s'1.  The  measured  response 
was  compared  with  predicted  ones  by  using  KHL  and  JC  models  (Fig.  3).  This 
experimental  response  was  not  included  in  the  determination  of  material  constants.  The 
predictions  were  made  assuming  that  10'5  s'1  span  is  under  isothermal  conditions  and  10'1 
s'1  and  1700  s'1  experiment  are  under  adiabatic  condition.  Especially  at  a  strain  rate  of 
1700  s'1,  and  initial  segment  at  10'5  s'1,  KHL  model  predictions  are  in  better  agreement 
with  the  measured  response  than  the  JC  model,  which  over  predicts  the  experimental 
results  at  these  two  strain  rates. 

4.2.2  Correlations  with  Nemat-Nasser  et  ai.  ’s  data  (Alloy  2) 

In  Fig.  4,  quasistatic  and  dynamic  measurements  (at  a  strain  rate  of  10'  s'  )  for 
various  temperatures  with  the  correlations  using  KHL  and  JC  models  are  shown.  These 
measurements  were  made  on  a  commercially  available  Ti-6A1-4V  alloy,  obtained  from 
Protech  Metals.  As  shown  in  Table  2,  the  composition  of  this  material  is  somewhat 
different  from  that  of  the  Ti-6A1-4V  alloy  studied  in  the  present  work  (Alloy  1),  see 
Table  1.  Since  the  quasistatic  data  to  determine  the  initial  material  constants  were  not 
available,  the  parameters  were  directly  optimized  using  the  parameters  obtained  for  the 
other  alloy  as  an  initial  input.  For  all  the  cases,  KHL  model  gives  better  correlation. 


4.2.3  Correlations  and  Predictions  Using  Macdougall  and  Harding ’s  data  (Alloy  3) 

The  data  from  torsional  experiments  were  extracted  and  converted  to  von  Mises 
effective  stress  and  effective  plastic  strain  assuming  a  =  and  £p  -y  i  4 3  where  r 
and  y  are  shear  stress  and  shear  strain,  respectively.  The  material  constants  determined 
from  the  correlation  with  three  torsional  measurements:  y  =  1000,  0.1,  and  0.0007  s'1 

(which  correspond  to  effective  von  Mises  strain  rates,  £  =  511  A,  0.0577,  and  0.0004  s'1, 
respectively).  These  constants  are  given  in  Tables  3  &  4.  Since  the  data  near  failure  may 
be  affected  by  necking  (tensile  loading),  or  buckling  in  case  of  torsional  loading,  these 
data  near  failure  are  not  included  in  the  analysis.  In  Fig.  5,  quasistatic  and  dynamic 
torsional  measurements  with  effective  strain  rates,  £  =  511  A,  0.0577,  and  0.0004  s'1  at  a 
temperature  of  293  K  are  shown  along  with  correlations  using  KHL  and  JC  models.  The 
material  constants  were  determined  from  the  same  experimental  results.  Again,  KHL 
model  correlations  are  better  than  those  with  JC  model  especially  at  higher  strain  rates. 
Based  on  the  material  constants  determined  only  from  the  torsion  experiments,  a 
comparison  between  tensile  dynamic  data  at  a  strain  rate,  £  =  2000  s'1,  and  predicted 
adiabatic  stress/strain  responses  using  KHL  and  JC  models  at  a  temperature  of  293  K  are 


shown  in  Fig.  6.  Again,  KHL  model  gives  better  prediction;  especially  the  thermally 
softened  work-hardening  due  to  heat  dissipation  is  more  closely  predicted  with  KHL 
model. 


Conclusions 


A  comprehensive  study  on  quasi-static  and  dynamic  responses  of  Ti-6A1-4V 
Titanium  (Alloy  1)  alloy  from  monotonous  and  strain-rate -jump  compression 
experiments  was  performed.  Quasistatic  experiments,  using  MTS  machine,  at  room 
temperature  were  conducted  on  cylindrical  titanium  specimens  at  strain  rates  of  10'  ,  10' 
and  1  s'1,  and  at  temperatures  of  233  K,  296  K,  422  K,  588  K,  and  755  K  at  a  strain  rate 
of  10'  s'  .  Dynamic  response  was  measured  using  the  split  Hopkinson  pressure  bar 
(SHPB)  technique  at  strain  rates  of  3378  s'1  at  room  temperature.  Experimental  results 
from  uniaxial  loading  experiments  were  used  to  determine  material  constants  for  the 
KHL  and  JC  constitutive  models  for  this  titanium  alloy.  It  was  found  that  Ti-6A1-4V  is 
nonlinearly  dependent  on  strain  rate,  as  well  as  temperature. 

The  KHL  model  was  found  to  correlate  better  than  JC  model  especially  during 
dynamic  deformation  regime.  The  thennal  softening  at  high  strain  rates,  together  with 
reduction  in  the  work  hardening  rate  with  increase  in  strain  rate  and  strain,  was  captured 
much  better  by  the  KHL  model  than  the  JC  model.  The  temperature  term  of  the  modified 
KHL  viscoplastic  constitutive  model  was  able  to  correlate  well  with  the  response  at  a 
temperature  (233  K)  which  was  lower  than  the  reference  temperature  (296  K);  the  JC 
model  was  not  valid  in  this  case. 

Material  constants  of  the  KHL  and  JC  models,  determined  previously  were  used  for 
predictions  of  the  observed  material  response  in  a  strain-rate-jump  experiment  to  examine 
the  validity  of  the  two  models.  The  test  was  performed  at  296  K  with  strain  rate  from  10'5 
s'1  to  10'1  s'1  and  then  to  1700  s'1.  Again  KHL  model  provided  better  agreement  with  the 
measured  response  than  the  JC  model. 

In  addition  to  this  alloy,  results  from  two  more  titanium  alloys  were  obtained  from  the 
existing  publications  (Nemat-Nasser  et  al.,  2001;  Macdougall  and  Harding,  1999)  and 
compared  with  correlations  and  predictions  by  using  KHL  and  JC  models.  Once  again, 
the  KHL  model  was  closer  to  the  observed  response  than  JC  model  especially  at  high 
strain  rates  and  at  different  temperatures. 
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Table  1:  Chemical  Composition  of  Ti-6AL-4V  (alloy  1)  in  the  present  investigation 


Element 

%  wt 

Composition 


Ti  Y 


Bal.  <0.0003 


N  C 


0.008  0.043 


H  Fe 


0.0041  0.15 


O  A1 


0.174  5.97 


V 


4.09 


Table  2:  Chemical  composition  of  commercial  Ti-6A1-4V  (alloy  2)  used  by  Nemat- 
Nasser  et  al.  (2001) 


Element  Ti  Y  N  C  H  Fe  O  Al  V 


^  .  Bal.  <0.001  0.01  0.01  0.0006  0.21  0.19  6.21  3.61 

Composition 


Table  3:  KHL  model  material  constants  determined  for  the  Ti-6A1-4V  (alloy  1)  and  the 
data  extracted  from  Nemat-Nasser  et  al.  (2001),  alloy  2,  and  Macdougall  and  Harding 
(1999),  alloy  3 


A 

(MPa) 

B 

(MPa) 

Hi 

n0 

C 

m 

Ti-6A1-4V  Alloy  1 

1069 

874.8 

0.5456 

0.4987 

0.02204 

1.3916 

Ti-6A1-4V  Alloy  2 

1097 

1004.7 

0.4993 

0.6268 

0.02219 

1.4796 

Ti-6A1-4V  Alloy  3 

1004 

325.1 

1.9015 

0.5263 

0.02204 

1.1206 

Table  4:  JC  model  material  constants  determined  for  the  Ti-6A1-4V  (alloy  1)  and  the  data 
extracted  from  Nemat-Nasser  et  al.  (2001),  alloy  2,  and  Macdougall  and  Harding  (1999), 
alloy  3 


A 

(MPa) 

B 

(MPa) 

n0 

C 

in 

Ti-6A1-4V  Alloy  1 
(with  optimization) 

1104 

1036 

0.6349 

0.01390 

0.7794 

Ti-6A1-4V  Alloy  1 
(without  optimization) 

1080 

1007 

0.5975 

0.01304 

0.7701 

Ti-6A1-4V  Alloy  3 

1119 

838.6 

0.4734 

0.01921 

0.6437 

Ti-6A1-4V  Alloy  4 

984 

520.3 

0.5102 

0.015 

0.8242 

Figure  1:  Quasi-static  and  dynamic  loading  experimental  results  (at  a  temperature  of 
296K)  for  different  strain  rates  with  correlations  using  KHL  and  JC  models;  the 
correlations  with  JC  model  were  shown  with  and  without  optimization  to  determine  the 
material  constants. 


3  1 

Figure  2:  Quasi-static  measurements  (at  a  strain  rate  of  10'  s'  )  for  various  temperatures 
with  correlations  using  KHL  and  JC  models;  the  correlations  with  JC  model  were  shown 
with  and  without  optimization  to  determine  the  material  constants.  Note  the  absence  of 
correlation  at  233  K  with  JC  model  in  which  (T-Tref)/(Tm-Tref)  becomes  negative. 
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Figure  3:  A  comparison  of  a  strain-rate -jump  experiments  results  with  predictions  using 
KHL  and  JC  models.  The  test  was  performed  at  296  K  with  strain  rate  from  10'5  s'1  to  10'1 
s'1  and  to  1700  s'1.  The  predictions  were  made  assuming  that  10‘5  s'1  region  is  under 
isothermal  conditions,  while  10'1  s'1  and  1700  s'1  regions  are  under  adiabatic  conditions. 
Material  constants  for  JC  model  used  here  were  determined  with  optimization. 
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Figure  4:  Quasi-static  (at  a  strain  rate  of  10'  s'  )  and  dynamic  measurements  for  various 
strain  rates  and  temperatures  (Nemat-Nasser  et  al.  2001)  with  correlations  using  KHL  and 
JC  model.  Material  constants  for  JC  model  used  here  were  determined  with  optimization. 
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Figure  5:  Quasi-static  and  dynamic  torsional  measurements  for  effective  strain  rates,  £  = 
511  A,  0.0577,  and  0.0004  s'1  at  a  temperature  of  293  K  (Macdougall  and  Harding,  1999) 
with  correlations  using  KHL  and  JC  models.  Material  constants  for  JC  model  used  here 
were  determined  with  optimization. 
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Figure  6:  Comparison  between  tensile  dynamic  data  with  a  strain  rate,  £'=2000  s'1,  and 
predicted  adiabatic  responses  using  KHL  and  JC  models  t  a  temperature  of  293  K.  The 
material  constants  were  optimally  determined  from  the  correlation  with  torsional 
experiments. 
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Abstract 

Results  from  a  series  of  experiments  on  three  different  titanium  alloys,  under 
quasi-static  and  dynamic  loading  conditions  are  presented.  The  Ti-6A1-4V  titanium  alloys 
include  the  ELI  version  and  two  with  higher  oxygen  content.  The  strain-rates  varied  from 
10‘6  s'1  to  3378  s'1  while  observations  were  made  at  temperatures  from  233K  to  755K. 
The  initial  and  deformed  microstructures  are  also  given.  Differences  in  these  alloys  are 
observed  in  terms  of  thermal  softening,  work  hardening,  and  strain-rate  &  temperature 
sensitivities.  The  Khan-Huang-Liang  (KHL)  model  is  used  to  effectively  simulate  the 
observed  response  obtained  from  these  experiments.  The  model,  with  constants 
determined  from  the  above  experiments,  is  then  used  to  predict  strain-rate  jump 
experimental  results  and  also  high  temperature  dynamic  experiments  for  one  of  the 
alloys;  the  predictions  are  found  to  be  very  close  to  the  observations. 
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Introduction 

Since  the  introduction  of  titanium  and  titanium  alloys  in  the  early  1950s  they  have 
become  important  materials  for  aerospace,  energy,  and  chemical  industries.  They  are 
used  chiefly  for  parts  that  require  good  corrosion  resistance,  moderate  strength  up  to  588 
K,  and  lightweight.  Titanium  alloys,  especially  Ti-6A1-4V,  an  a+(3  type  titanium  alloy  is 
a  largely  used  alloy  in  many  industries  because  of  its  extremely  attractive  properties  like 
high  specific  strength,  good  deformability,  reasonable  ductility  and  ability  to  withstand 
high  temperatures  and  resistance  to  corrosion.  It  is  primarily  used  in  aero-engine,  gas 
turbines  and  other  applications.  The  development  of  a  relatively  economical  Ti-6A1-4V 
alloy,  with  a  low  interstitial  content  has  prompted  a  lot  of  interest  in  its  possible  use  in 
annor  tanks  because  of  improved  ductility,  whereas  the  conventional  more  expensive  Ti- 
6A1-4V  alloy  has  been  used  primarily  in  aerospace  components  chiefly  because  of  its 
high  strength  to  weight  ratio.  This  advantage  has  also  been  contemplated  in  various 


applications  like  armor,  including  ceramic  tiles  encapsulated  in  these  titanium  alloys. 
Several  experimental  studies  have  been  carried  out  on  similar  materials,  which  included 
quasi-static  and  dynamic  experimentation  (Follansbee  and  Gray,  1989;  Lee  and  Lin,  1997 
&  1998;  da  Silva  and  Ramesh,  1997;  Lesuer,  2000;  Nemat-Nasser  et  al.,  2001;  Majorell 
et  al.,  2002). 

Since  the  mechanical  properties  of  this  alloy  are  heavily  influenced  by  its  impurities 
and  heat  treatment  it  is  imperative  that  this  alloy  be  understood  in  terms  of  the 
microstructure,  constituents,  and  the  processing  history.  The  unalloyed  titanium  exists  in 
mainly  two  forms,  a-titanium  (hep  structure)  at  room  temperature,  and  at  high 
temperatures  it  exists  as  (3-titanium  (bee  structure).  The  a  phase  (primary  and  secondary) 
is  more  dominant  in  the  temperature  ranging  from  300  K  to  800  K,  while  the  amount  of  P 
phase  starts  to  increase  thereon  and  is  almost  totally  P-phase  at  about  1270  K,  which  is 
called  the  P  transus  temperature  [Majorell  et  al.,  2002],  The  a  phase  is  much  harder  than 
the  P  phase,  so  the  drop  in  stress  levels  upon  increase  of  temperatures  can  be  noticed  with 
increasing  P  content.  However,  formability  increases  as  the  phase  changes  from  a  to  p. 
Aluminum  and  Vanadium  are  added  to  titanium  to  stabilize  the  two  phases.  Aluminum, 
oxygen,  nitrogen,  and  carbon  are  the  a  stabilizers,  while  vanadium,  iron  and  manganese 
are  the  P  stabilizers.  Planar  slip  is  more  dominant  in  the  lower  temperature  range  (<500 
K).  The  slip  occurs  on  the  basal,  prismatic  and  the  pyramidal  planes.  At  temperatures 
below  800  K,  AI  transforms  the  primary  a  into  secondary  a  precipitates.  The  secondary  a 
phase  has  a  higher  dislocation  density  than  primary  a  and  has  a  smaller  grain  size.  This 
basically  suggests  that  it  is  harder  than  primary  a  phase  [Picu  et  al. ,  2002],  Also,  studies 
by  Conrad  and  Wang  (1979),  where  they  studied  high  purity  titanium  against  titanium 
with  7.4%  Al,  suggested  an  increase  in  the  flow  stress  with  the  addition  of  Al  content. 
This  is  because  the  increase  in  the  presence  of  Al  concentration  enhances  the  production 
of  secondary  a  phase,  which  consequently  makes  the  alloy  harder.  They  also  showed  that 
there  was  an  inverse  dependency  of  stress  to  the  grain  size  for  the  same  material.  The 
effect  of  grain  size  on  the  yield  and  flow  stress  of  titanium  was  found  to  obey  the  Hall- 
Petch  relation  (Conrad,  1984).  Also,  the  strain  hardening  of  the  alloy  is  more  due  to  the 
presence  of  the  lamellar  Widmanstatten  microstructure  of  the  secondary  a  phase. 
However,  Sargent  and  Conrad  (1970)  noticed  that  the  addition  of  aluminum  does  not  alter 
the  rate  controlling  mechanism  and  the  aluminum  contributes  only  to  the  athermal 
components  of  the  flow  stress  whereas  the  rate  controlling  mechanism  is  the  thermally 
activated  component  overcoming  of  the  interstitial  solute  atoms.  De  Meester  and  Doner 
and  Conrad  (1975),  found  that  between  temperatures  from  300  K  to  800  K  the  flow  stress 
vs.  temperature  curve  of  a-titanium  was  almost  parallel  to  that  of  Ti-6A1-4V  alloy  of 
nearly  same  interstitial  content,  concluding  that  alloying  elements  had  nothing  to  do  with 
temperature  sensitivity  in  this  range  and  that  the  interaction  of  dislocations  with 
interstitial  impurities  (C,  N,  O  and  H)  were  more  responsible.  They  also  concluded  that 
the  deformation  kinetics  and  the  related  parameters  for  the  Ti-6A1-4V  alloy  are  the  same 
as  those  for  unalloyed  titanium. 

Chichilli  et  al.  1998  found  that  twinning,  an  important  dislocation  mechanism, 
occurred  at  all  temperatures  below  773  K  for  a-titanium;  at  high  strain  rates,  the  density 
of  twins  significantly  increased.  There  was  no  twinning  seen  for  Ti-6A1-4V  alloys  even  at 
low  temperatures  of  100  K  (Paton  et  al.  1976),  suggesting  that  increase  in  aluminum 


content  inhibited  twinning.  However,  at  room  temperature,  there  was  some  twinning 
observed  by  Follansbee  and  Gray  (1989)  at  strain  rate  of  5000  s'1  in  grains  which  were 
larger  than  the  average  size. 

Conrad  et  al.  (1975)  proposed  an  equivalent  oxygen  content  (Oeq  =  O+2N+0.75C  ), 

which  gives  the  effect  of  dislocations-impurity  interaction  on  the  yield  strength  of  the 
material.  They  also  noted  that  the  strengthening  of  the  alloy  due  to  these  interstitials  was 
dependent  on,  in  decreasing  order,  Carbon,  Oxygen  and  followed  by  Nitrogen.  Effect  of 
Hydrogen  in  most  cases  was  neglected.  Majorell  et  al.  (2002)  studied  the  deformation 
characteristics  and  response  of  a  textured  and  untextured  Ti-6A1-4V  alloy  rod  that  was 
manufactured  by  Allvac.  This  study  was  over  a  strain-rate  range  of  10'  to  10  s'  and  a 
temperature  range  of  650-1345  K  (710-1970  F).  They  did  not  observe  any  dynamic 
strain  aging  at  any  of  the  temperatures  or  strain  rates  investigated.  The  equivalent  oxygen 
content  of  the  specimen  was  0.206%  and  the  average  grain  size  of  the  specimen  was  20 
pm.  They  also  included  the  equivalent  oxygen  content  and  aluminum  impurities  in  the 
thermal  stress  component  and  the  grain  size  effect  in  the  athermal  stress  component. 

The  aim  of  this  study  is  to  compare  the  responses  of  the  three  alloys  and  subsequently 
present  the  modeling  capability  of  the  KHL  model  (Khan  and  Liang,  1999),  which  was 
modified  further  in  an  earlier  paper  (Khan  et  al.  2004),  for  three  titanium  alloys. 
Comparisons  of  KHL  model  correlation  and  predictions  with  Johnson-Cook  model 
(Johnson  and  Cook,  1983)  were  also  given  in  that  paper  for  one  of  the  alloys  presented 
here,  and  also  alloys  used  by  Nemat-Nasser  et  al.  (2001),  &  Macdougall  and  Harding 
(1999).  The  reader  is  referred  to  papers  by  Abed  and  Voyiadjis  (2005)  &  Uenishi  and 
Teodosiu  (2004),  for  other  modeling  approaches. 


Experimental  Procedure 

2. 1  Material 

Three  titanium  alloys  with  different  compositions  are  included  in  this  study  of  the 
effect  of  impurities  on  the  stress-strain  response  of  the  material  under  different  loading 
conditions.  The  alloys  were  in  as  received  condition  in  the  form  of  plates,  and  specimens 
were  machined  to  a  diameter  of  0.5  in.  (12.7mm)  and  length  0.75  inches  (19.05  mm).  The 
dynamic  specimens  were  machined  to  discs  with  0.4  in.  (10.16  mm)  diameter  and  0.2  in 
(5.08  mm)  thickness.  The  loading  direction  in  all  these  specimens  was  the  thickness 
direction  (axis  of  cylinders)  of  the  plate.  The  chemical  composition  of  these  three 
different  alloys  is  given  in  Table  1.  These  alloys  have  been  designated  with  different 
numbers  i.e.  as  alloys  1,  2  and  3.  The  equivalent  oxygen  content  for  alloys  1,  2  and  3 
were  found  to  be  0.229,  0.152  and  0.222  respectively.  Alloys  1  &  3  were  manufactured 
using  single  electron  beam  melting  process,  resulting  in  higher  oxygen  content  than  the 
ELI  version  (alloy  2). 


2.2  Microstructure  Evaluation 


The  microstructure  of  these  alloys  was  determined  before  and  after  deformation  using 
a  scanning  electron  microscope  (SEM).  The  specimens  were  cut  along  the  compression 
axis  of  the  specimen  (the  thickness  direction  of  the  plate).  The  preparation  of  the 
specimens  involved  polishing  using  a  Buehler  Metaserv  grinding-polishing  machine. 
Different  grits  of  carbide  paper  were  used  for  grounding  starting  from  240,  400,  600,  800, 
to  2400.  The  specimens  were  then  polished  again  using  a  Struers  DP-Plus  micro  cloth 
with  slurry  of  0.3  pm  alumina  to  get  a  mirror-like  finish.  It  was  eventually  etched  with 
KrolTs  reagent  (a  solution  of  1%  HF,  2%  HNO3  and  97%  distilled  HoO),  for  about  5  min. 

The  photomicrographs  of  the  initial  and  deformed  microstructures  revealed  an  alpha- 
beta  type  microstructure,  which  is  typical  for  the  alloy  (Fig.  1  and  2).  It  confirms  that  the 
alloys  went  through  different  annealing  and  cooling  processes  with  higher  amount  of 
transfonned  P  content  in  alloy  3.  This  alloy  has  gone  through  annealing  at  higher 
temperatures  (~  1750F).  Its  microstructure  consists  of  the  hep  a  grains  and  a  transformed 
P  with  a  lamellar  Widmanstatten  structure  dispersed  between  these  primary  globular  a 
grains.  The  higher  annealing  temperature  along  with  higher  rate  of  cooling,  is  responsible 
for  the  increased  transformed  P  content  as  compared  to  the  other  alloys.  The  percentage 
of  lamellar  transfonned  P  is  lower  in  alloy  2  relative  to  alloy  3.  However,  the  aluminum 
content,  which  is  an  a  phase  stabilizer,  is  higher  in  alloys  1  and  2  than  alloy  3.  Also,  the 
average  grain  size  for  the  titanium  alloys  was  approximately  9  pm  as  calculated  using 
intercept  method  (ASTM  El  12-96). 

The  deformed  photomicrographs  of  alloys  2  and  3  are  shown  in  Figures  1(b)  and  2(b), 
respectively.  The  specimen  chosen  for  the  microstructural  examination  were  deformed  at 
somewhat  higher  temperatures  (422  K)  for  alloys  3  than  in  case  of  alloy  2,  which  was 
deformed  at  room  temperature  (296  K).  The  strain  rate  during  the  deformation  of  alloy  2 
was  1  s'  whereas  for  alloy  3  it  was  10'  s'  .  The  deformed  micro  structures  in  both  cases 
show  elongated  hep  a  grains.  The  P  grains  in  the  samples  displayed  no  apparent  change 
in  the  substructure  and  this  was  same  for  all  the  alloys  independent  of  the  level  of 
deformation  or  loading  conditions.  This  was  consistent  with  the  other  researchers 
(Follansbee  and  Gray,  1989).  Also,  very  little  to  almost  no  twinning  was  observed  in  the 
alloys.  This  is  also  consistent  with  other  studies  (Follansbee  and  Gray,  1989),  where  the 
authors  found  deformation  twins  only  at  very  high  strain  rate  (5000  s'1)  experiments.  The 
reason  behind  this  behavior  is  that  at  very  high  strain  rate  or  at  very  low  temperatures, 
higher  stresses  are  required  to  cause  plastic  deformation  and  since  the  dislocation 
mobility  is  low  at  these  conditions,  the  deformation  twins  develop  as  a  consequence. 
However,  in  the  case  of  low  strain  rate  and  high  temperatures,  dislocation  mechanisms 
are  supposed  to  be  the  dominating  components  in  the  deformation  mechanism  of  the 
alloys  through  planar  slips  with  almost  negligible  amounts  of  deformation  twins. 

The  photomicrographs  of  deformed  samples  under  dynamic  loading  conditions  at 
various  temperatures  for  alloy  3  are  also  shown  in  the  Fig  3.  The  specimens  are  cut  along 
the  compression  axis  for  this  investigation.  The  samples  show  that  elongation  of  grains 
occurs  for  the  deformed  sample  perpendicular  to  the  loading  axis.  The  grain  elongation 


ratio,  as  per  ASTM  standard  El  12-96,  for  this  alloy  was  calculated  for  undeformed  and 
deformed  samples  and  is  shown  in  Table  3.  Also,  the  average  grain  size  for  each  of  the 
photomicrographs  for  this  alloy  are  calculated  and  shown  in  the  table.  The  average  grain 
size  was  found  to  decrease  with  increasing  strain  rates.  It  is  believed  that  for  high  strain 
rate  and  high  temperature  deformations,  the  energy  generated  causes  dynamic 
recrystallization  of  the  platelet  type  alpha  into  globular  primary  alpha  grains  (422K  and 
715K)  as  the  platelet  type  alpha  content  is  lowest  in  these  samples. 

2.3  Quasi-Static  Experiments 

2.3.1  Room  temperature  experiments 

The  quasi-static  experiments  perfonned  were  in  compression.  Cylindrical  specimens 
were  made  out  of  the  base  plates  into  dimensions  of  approximately  0.75  in.  (19.05  mm) 
in  length  and  0.5  in.  (12.7  mm)  diameter  with  the  longitudinal  axes  in  the  thickness 
direction.  Care  was  taken  during  the  machining  of  the  specimens  that  no  plastic  work 
hardening  is  introduced  on  the  surface  and  the  last  three  machining  cuts  were  one 
thousands  of  an  inch.  These  specimens  were  then  subjected  to  different  strain  rates 
ranging  from  10'6  to  1  s'1  on  the  MTS  servo  hydraulic  axial/torsional  material  testing 
machine.  The  1  s'1  experiment  was  performed  as  a  load-unload-reload  experiment  with 
increments  of  5%  strains  with  a  30-minute  interval  between  each  consecutive  loading,  to 
obtain  an  isothermal  response  of  the  material.  For  these  room  temperature  compression 
experiments,  KFEL-2-120-C1  high  elongation  uniaxial  strain  gages,  manufactured  by 
Kyowa  Ltd.  (Japan),  were  mounted  diametrically  opposite  to  each  other  at  the  center  of 
each  specimen  using  an  AE-15  (Micro  Measurements  Inc.)  epoxy  resin.  The  gages  were 
then  connected  to  the  MTS  TestStar  software  through  signal  conditioners.  The  strains 
obtained  using  strain  gages  were  corrected  before  conversion  into  true  plastic  strain  for 
Wheatstone  bridge  non-linearity  and  gage  factor  variation  (Huang  and  Khan,  1991).  The 
MTS  transducers  directly  supplied  load  data  to  the  software.  The  interface  between  the 
fixtures  and  the  test  specimens  for  each  of  the  compression  experiments  was  lubricated 
with  Teflon  sheets  and  a  high  vacuum  grease  in  order  to  achieve  uniformity  in 
deformation  and  avoid  barreling  of  the  specimen,  maintaining  a  uniform,  uniaxial  stress 
state. 

2.3.2  Experiments  at  different  temperatures 

The  stress-strain  response  was  observed  during  experiments  at  temperatures  of  233, 
422,  589,  and  755  K  at  a  strain  rate  of  10‘4  s'1  for  alloys  1  and  2,  whereas  for  alloy  3, 
strain  rate  of  10"  s'  was  used  to  perform  these  experiments.  Alumina  bars  were  attached 
to  the  grips  of  the  MTS  using  Vascomax  C-350  steel  fixtures.  The  specimen  was  placed 
between  the  alumina  bars.  A  thermocouple  was  cemented  to  the  center  of  the  specimen  to 
monitor  the  temperature  throughout  the  experiment.  Once  the  desired  temperature  was 
reached,  it  was  kept  constant  for  30  minutes  so  that  steady  state  conditions  existed  and 
consequently  ensured  a  uniform  temperature  throughout  the  specimen.  Extreme  high 
temperature  grease  was  applied  between  the  specimen  and  the  ceramic  bars  to  reduce 
frictional  effects  and  have  uniaxial  stress  state  throughout  the  length  of  the  experiment. 


The  LVDT  of  the  MTS  supplied  stroke  data  in  the  fonn  of  displacement  and  force  was 
supplied  by  the  force  transducers.  The  displacement  data  supplied  by  the  machine 
included  machine  elastic  deformations,  the  deformation  of  the  ceramic  bars  and  the 
variation  in  thickness  of  the  layer  of  grease.  This  measurement  in  displacement  was 
corrected  by  running  a  “blank  experiment”,  i.e.  without  the  specimen  and  keeping  all  the 
other  parameters  approximately  same.  The  displacement  for  the  blank  experiment  was 
then  subtracted  from  the  experiment  displacement  at  each  load  to  get  the  actual 
displacement  of  the  specimen.  The  remaining  calculations  were  done  using  this  corrected 
displacement. 

For  the  low  temperature  experiments,  dry  ice  in  form  of  small  granules  was  used.  A 
chamber  was  placed  around  the  specimen  with  dry  ice.  The  temperature  of  the  specimen 
was  monitored  using  the  same  technique  as  in  the  case  of  high  temperature  experiments. 
The  specimen  temperature  went  down  to  233  K  and  remained  steady  during  the  test.  Low 
temperature  grease  was  used  as  lubricant  between  the  specimen  and  the  alumina  bars. 
Corrections  for  the  displacement  were  performed  similar  to  that  of  high  temperature 
experiments 

2.4  Dynamic  experiments  at  room  temperature 

Dynamic  experiments  on  the  alloy  were  performed  using  the  compression  split- 
Fiopkinson  pressure  bar  technique  (SHPB),  which  is  a  widely  used  technique  for  high 
strain-rate  experiments.  Details  of  this  technique  can  be  found  in  Khan  and  Liang  (1999). 
The  diameter  of  these  bars,  made  from  Vascomax  C-350  steel,  was  a  constant  0.5  in. 
(12.7  mm).  The  specimen  was  machined  to  the  shape  of  a  disc  and  was  sandwiched 
between  the  incident  and  transmitted  bars.  The  dimensions  of  the  specimens  used  for 
dynamic  experiments  were  0.2  in.  (5.1  mm)  in  thickness  and  0.4  in.  (10.2  mm)  in 
diameter.  The  interfaces  of  these  incident  and  transmitted  bars  with  the  specimen  were 
lubricated  with  grease  to  reduce  friction  and  to  maintain  a  uniform  stress  state  in  the 
specimen.  The  strain  gages  mounted  on  the  bars  were  manufactured  by  Kyowa  Inc., 
Japan.  These  gages  were  bonded  diametrically  opposite  to  each  other  using  an  epoxy 
resin.  The  position  of  the  gages  was  such  that  there  was  no  interaction  of  the  incident, 
reflected  and  transmitted  waves  with  each  other.  These  gages  were  connected  to  a  Nicolet 
440  oscilloscope  through  a  potentiometer  circuit.  The  stresses  and  the  corresponding 
strains  in  the  specimen  are  calculated  using  the  incident,  reflected  and  transmitted  waves 
signals.  More  detailed  information  about  the  split-Hopkinson  can  be  found  elsewhere 
(Follansbee,  1979). 

For  high  temperature  deformations,  the  disk  specimens  were  heated  to  a 
predetennined  temperature.  Thermocouples  were  cemented  to  the  outside  surface  of  the 
specimens  to  monitor  the  current  temperatures.  Once  the  desired  temperature  was 
achieved,  it  was  held  constant  for  about  15  minutes  at  that  value  to  insure  uniformity  of 
heating  at  the  surface  and  the  core  of  the  samples.  The  Hopkinson  bars  were  kept  cooler 
using  cooling  devices.  The  temperature  in  the  bars  was  kept  constant  around  room 
temperature  so  as  not  to  cause  any  significant  change  in  the  elastic  modulus. 


During  high  strain  rate  experiments  most  of  the  heat  generated  by  the  plastic  work 
is  entrapped  within  the  specimen  making  it  somewhat  relatively  soft,  this  process  is 
called  thermal  softening.  Thermal  softening  is  a  function  of  the  specific  heat  capacity  at 
constant  pressure,  mass  density  and  the  deformation  in  the  specimen.  The  rest  of  the  heat 
generated  is  used  up  to  change  the  microstructure,  e.g.  increase  the  dislocation  density. 
Thus,  high  strain-rate  deformation  is  adiabatic  and  there  is  a  significant  increase  in  the 
temperature  of  the  specimen  during  deformation.  The  effect  of  this  softening  is 
significant  and  cannot  be  neglected.  The  temperature  rise  in  the  specimen  can  be 
calculated  by  the  following  equation: 

AT  =  — —  frxleWf  (1) 
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In  this  study,  the  fraction  of  plastic  work  converted  to  heat  is  taken  to  be  0.9,  which 
implies  that  90%  of  the  plastic  work  done  during  the  deformation  is  assumed  to  be 
converted  to  heat  and  is  responsible  for  the  rise  in  the  temperature  of  the  specimen,  p  and 
Cp  are  the  mass  density  and  the  specific  heat  at  constant  pressure  for  the  alloy.  This 
relation  is  used  to  convert  the  original  experimental  response  at  high  strain  rates  under 
adiabatic  heating  to  corresponding  isothermal  behavior.  The  density  for  the  titanium  (Ti- 
6A1-4V)  alloys  is  4428  kg/m3  (Lesuer,  2000)  and  the  specific  heat  capacity  at  constant 
pressure  for  the  alloys  can  be  expressed  as  a  function  of  the  temperature  of  the  material 
(Military  Handbook,  1998). 

C=  559.77  -  0.1473T  +  0.00042949T2  J/(kg-K)  (278  K  <  T  <  1 144  K).  (2) 


2.5  Compressive  Strain  rate  jump  experiments 

Compressive  strain  rate  jump  experiments  were  performed  on  the  titanium  alloys 
as  a  part  of  this  comprehensive  study.  The  quasi-static  regime  part  of  the  experiment  was 
initially  perfonned  on  the  MTS.  The  strain  rate  of  the  samples,  once  loaded  to  a  certain 
level  of  strain,  was  changed  abruptly  and  the  corresponding  response  of  the  alloy  was 
recorded.  After  the  predetermined  level  of  strain  was  achieved,  the  sample  was  unloaded 
and  machined  to  a  disc  form  for  a  dynamic  compression  experiment  using  the  split- 
Hopkinson  pressure  bar  technique. 


2.  Modeling  Procedure 

The  KHL  model  utilizes  a  systematic  method  of  determining  initial  material  constants 
(see  Khan  and  Liang,  1999),  which  are  used  as  input  in  a  program  to  obtain  refined 
material  constants  by  employing  an  optimization  technique  and  to  predict  the  material 
response  over  a  wide  range  of  strain  rates  and  temperatures.  There  are  six  different 
material  constants,  which  define  the  material  behavior  with  respect  to  yield  strength, 
strain  hardening,  strain  rate  sensitivity,  temperature  sensitivity  and  strain  hardening 
sensitivity  to  strain  rate.  The  model  is  more  suited  for  modeling  the  material  response  as 
compared  to  JC  model,  which  fails  to  provide  observed  decreasing  work  hardening 


behavior  in  some  materials  with  increasing  strain  rate.  More  detailed  advantages  of  this 
model  over  JC  model  are  given  in  Khan  et  al.  (2004).  The  KHL  constitutive  equation  is 
given  as: 
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where,  a  is  the  true  (Cauchy)  stress  and  ep  is  the  true  plastic  strain.  Tm,  T ,  Tref  are 

melting,  current,  and  reference  temperatures,  respectively.  D£=  106  s'1  (arbitrarily  chosen 

upper  bound  strain  rate)  and  £=\  s"1  (reference  strain  rate,  at  a  reference  temperature  of 
Tref,  usually  room  temperature,  at  which  material  constants  A,  B  and  no  are  determined). 
e  is  the  current  strain  rate,  n , ,  C  and  m  are  additional  material  constants.  For  Ti-6A1-4V 
alloys,  the  melting  temperature  was  taken  to  be  1933  K  (Nemat-Nasser  et  al.,  2001).  The 
reference  temperature  was  the  constant  room  temperature  for  experiments  at  296K. 


The  material  constants  for  this  model  are  determined  using  several  uniaxial  quasi¬ 
static  and  dynamic  room  temperature  experimental  results  at  different  strain  rates  and 
high  &  low  temperature  experimental  results  at  one  strain  rate.  This  usually  provides  one 
with  a  set  of  initial  material  constants  which  can  be  used  as  input  into  a  programming 
software  which  utilizes  constrained  optimization  procedure  using  least  square  method  to 
obtain  a  set  of  more  refined  material  constants.  The  dynamic  data  before  being  input  in 
the  optimization  scheme  is  converted  into  an  equivalent  isothermal  data  by  calculating 
the  instantaneous  temperature  rise  in  the  material  at  any  strain  level  using  Eq  (1)  and 
calculating  stress  increment  for  that  temperature  rise  to  get  corresponding  isothermal 
stress  response.  The  material  constants  thus  determined  were  used  to  predict  a  strain  rate 
jump  experiment,  which  was  performed  separately  for  each  alloy.  This  experiment  was 
not  included  in  the  calculation  of  the  material  constants.  The  strain  rate  jump  experiment 
included  at  least  two  different  strain  rates  in  the  quasi-static  regime  followed  by  a 
dynamic  regime  strain  rate. 


3.  Results  and  correlations  with  the  constitutive  model 


The  comprehensive  observed  response  of  the  alloys  are  given  by  symbols  (e.g. 
v,0,x,+,o,  etc.),  in  Fig  4-13.  For  alloy  1  five  different  strain-rate  experiments  were 
performed  ranging  from  10~6  s'1  to  dynamic  strain  rate  of  3100  s'1.  The  material  behavior 
was  found  to  be  strain  rate  sensitive  and  has  a  non  linear  work  hardening  response  at  low 
deformation  levels  transitioning  to  an  almost  linear  response  at  higher  deformation  levels. 
Also,  there  is  a  slight  reduction  in  work  hardening  rate  with  increasing  strain  rates  in  the 
quasi-static  loading  regime.  The  JC  model,  with  no  provisions  to  replicate  such  behavior 
was  found  earlier  not  as  good  as  the  corresponding  KHF  model  (Khan  et  al.,  2004).  There 
is  no  or  little  work  hardening  effect  seen  in  the  adiabatic  dynamic  experiment  due  to 
thermal  softening.  The  observed  response  of  the  material  is  also  a  non-linear  function  of 
temperature.  The  increase  in  test  temperature  causes  a  drop  in  the  flow  stress.  A  slight 
decreasing  work  hardening  rate  behavior  can  be  noticed  as  the  temperature  increases  and 
also  with  an  increase  in  the  level  of  strain.  This  is  consistent  with  other  studies  on  the 
alloy  (Majorell  et  al.,  2002,  Fee  and  Fin,  1997,  etc.).  This  kind  of  behavior  is  amplified  in 


the  experiment  performed  at  755  K  temperature.  The  flow  stresses  in  the  alloy  is  also 
found  to  be  more  sensitive  to  temperature  of  the  material  than  to  the  strain  rate  as  also 
found  by  Nemat-Nasser  et  al.  (2001)  in  a  similar  alloy.  For  alloy  2  experiments  were 
performed  at  strain  rates  ranging  from  10'6  s'1  to  2926  s'1.  A  low  temperature  experiment 
at  233  K  was  also  performed  for  this  alloy.  Other  quasi-static  experiments  at  different 
temperatures  were  performed  at  422  K,  588  K  and  755  K.  All  of  these  experiments  were 
well  below  beta  transus  temperature  (-1270K)  to  ensure  that  the  material  did  not  change 
phases  during  the  experiments.  Similarly,  different  strain-rate  and  temperature 
experiments  were  performed  for  alloy  3.  For  room  temperature  experiments  the  range  of 
strain  rate  was  from  10'5  s'1  to  3378  s'1.  Experiments  at  different  temperatures  (233K, 
422K  588K  and  755K)  on  this  alloy  were  perfonned  at  a  strain  rate  of  10'  s'  . 


4.1  Correlations  and  predictions  for  alloy  1 

The  correlations  obtained  for  room  temperature  experiments  are  also  shown  in  Fig.  4 

using  lines  ( - ,  — ,  -  etc.  ).  The  various  material  constants  obtained  for  this  alloy 

using  MATLAB  software,  utilizing  the  least  square  method  with  constrained 
optimization,  are  shown  in  Table  2.  The  constants,  depicting  the  strain  rate  sensitivity  of 
work  hardening  nj  and  the  strain  rate  sensitivity  parameter  C,  were  kept  same  for  all  three 
alloys.  In  spite  of  putting  these  extra  constraints  in  the  optimization  procedure  the  KHL 
model  was  able  to  successfully  correlate  with  the  material  response  over  large  strains  and 
displaying  excellent  capability  in  simulating  decreased  work  hardening  behavior  with 
increasing  strain  rate.  The  various  correlations  for  different  temperatures,  using  lines,  can 
be  seen  in  Fig.  5.  Here  too,  the  model  correlates  reasonably  well  with  the  data.  Further,  in 
Fig.  6,  a  strain  rate  jump  experiment  from  10'5  s'1  to  10"1  s'1,  followed  by  another  jump  at 
a  dynamic  strain  rate  of  1816  s'1,  is  predicted  by  the  model.  This  experiment  was 
performed  at  room  temperature.  These  results  were  not  used  in  the  calculation  of  the 
material  constants,  earlier.  The  predictions  from  the  model  are  shown  by  a  solid  line;  the 
model  predicts  the  observed  response  very  closely  in  this  strain  rate  jump  experiment. 

4.2  Correlations  and  predictions  for  alloy  2 

The  various  material  constants  for  this  alloy  are  given  in  Table  2.  Here  again,  the 
two  material  constants,  «/  and  C,  were  kept  same  as  it  was  assumed  that  the  heat 
treatment  and  slight  changes  in  composition  only  changes  constants  A,  B,  no  and  m. 
Results  of  the  correlations  can  be  seen  in  Fig.  7.  The  model  is  in  good  agreement  with  the 
experimentally  observed  response  during  the  quasi-static  room  temperature  experiments 
throughout  the  deformation  regime.  The  dynamic  response  is  also  fairly  close  to  the 
correlations,  including  the  replication  of  the  reduced  work  hardening  throughout  the 
deformation.  The  model  correlations  for  observed  responses  at  different  temperatures  are 
given  in  Fig.  8.  The  model  is  again  in  good  agreement  with  the  experimentally  observed 
responses  for  the  low  (233  K)  as  well  as  other  high  temperature  experiments.  Similar  to 
alloy  1,  a  strain  rate  jump  experiment  was  performed  on  this  alloy  from  10'5  s'1  to  1  s'1 
and  then  back  to  10'5  s'1,  before  performing  a  dynamic  experiment  on  the  same  specimen. 


Again,  this  experiment  was  not  used  to  determine  the  material  constants.  The  model 
predictions  are  in  good  agreement  for  all  strain  rates  and  to  all  levels  of  strain  as  shown  in 
Fig  9.  In  predictions,  the  responses  were  assumed  adiabatic  for  the  strain  rate  of  1  s'1  and 
also  the  dynamic  strain  rate  regime  portion  of  the  experiment. 

4.3  Correlations  and  predictions  for  alloy  3 

The  response  modeling  for  alloy  3  was  perfonned  similar  to  the  other  two  alloys.  The 
material  constants  thus  obtained  after  refinement  of  the  initial  constants  are  given  in 
Table  2.  The  measured  responses  during  four  room  temperature  quasi-static  experiments, 
and  the  dynamic  experiment  also  at  room  temperature  are  shown  in  Fig.  10  along  with  the 
correlations  obtained  for  this  alloy.  The  correlations  are  in  excellent  agreement  with  the 
observed  responses.  For  the  dynamic  case,  the  isothermal  values  were  converted  into 
adiabatic  response  before  plotting  for  comparison  with  the  adiabatic  dynamic  results. 
This  was  done  in  the  cases  of  other  two  alloys  as  well.  The  correlations  with  the 
experimental  results  at  different  temperatures  and  at  the  strain  rate  of  10'  s'  are  provided 
in  Fig.  11.  The  observations  and  the  model  correlations  are  again  in  good  agreement  with 
each  other  at  all  temperatures.  Furthermore,  as  in  the  case  of  the  other  two  alloys,  a  strain 
rate  jump  experiment  was  perfonned  at  a  strain  rate  of  10'5  s'1  to  10'1  s'1  followed  by  a 
dynamic  strain  rate  experiment  at  1700  s'1  on  the  same  specimen.  Predictions  were 
performed  based  on  the  material  constants  obtained  for  the  alloy  from  other  experimental 
results.  The  experimental  results  along  with  the  predictions  are  compared  in  Fig.  12.  As, 
in  all  the  previous  cases  the  predictions  are  in  good  agreement  with  the  observed  results. 
The  predictions  for  the  strain  rate  of  10'1  s'1  were  converted  into  corresponding  adiabatic 
response,  before  comparison,  as  there  was  a  temperature  increase  noticed  in  the  specimen 
during  the  experiment.  Dynamic  strain  rate  regime  predictions  are  also  converted 
similarly  in  order  to  obtain  the  adiabatic  response. 

Predictions  for  high  strain  rate  experiments  are  also  shown  for  different 
temperatures  in  Fig.  13.  These  predictions  were  also  converted  to  adiabatic  response  and 
plotted  along  with  the  experimental  results.  The  predicted  response  is  again  in  excellent 
agreement  with  the  experimental  observations. 

Overall,  the  correlations  and  predictions  using  are  very  close  to  measured 
experimental  responses  for  a  large  range  of  strain  rates  and  temperatures.  Further,  it  was 
observed  from  the  measured  responses  that  alloy  1  and  3  had  higher  flow  stress  as 
compared  to  the  alloy  2.  This  is  attributed  to  the  fact  that  the  interstitial  solute  content  or 
the  equivalent  oxygen  content  (Oeq),  which  is  known  to  strengthen  the  alloy,  is  higher  in 
these  alloys  (Conrad  et  ah,  1975).  Alloy  2  is  the  ELI  grade  alloy  and  has  a  lower  flow 
stress  response  during  the  various  experiments.  The  work  hardening  in  the  alloys  1  and  3 
were  found  to  be  quite  similar  but  higher  than  that  of  alloy  2  and  this  is  also  confirmed  in 
Table  2,  where  the  material  constants  ( B  and  no  )  are  listed.  However  there  is  a  slightly 
different  response  of  the  alloys  at  different  temperatures,  this  characteristic  is  also 
highlighted  in  Table  2,  where  the  three  alloys  have  different  temperature  sensitivity 
parameter.  From  the  above  observations  it  can  be  concluded  that  the  interstitial  impurities 
have  significant  effect  on  the  flow  stress  and  work  hardening  response  of  the  alloys,  and 


these  additions  also  affect  the  temperature  sensitivity  of  the  alloys.  This  is  consistent  with 
the  findings  of  other  researchers  (de  Meester  and  Doner  and  Conrad,  1975).  Also,  the 
interstitial  content  has  no  apparent  effect  on  the  work  hardening  sensitivity  of  strain  rate, 
as  well  as  the  overall  strain  rate  sensitivity  of  the  material. 

Conclusions 

In  the  present  paper,  the  quasi-static  and  dynamic  responses  of  the  three  titanium 
alloys  were  studied  in  a  systematical  manner  over  large  ranges  of  strain  rates  and 
temperatures.  The  alloys  were  found  to  be  more  sensitive  to  change  in  temperature  than 
strain  rate.  The  KHL  model  was  utilized  to  model  the  observed  behavior  of  these  three 
different  alloys  over  this  range.  The  model  was  found  to  be  fairly  accurate  to  capture 
closely  all  the  important  features  observed  for  these  alloys.  Material  constants  of  the 
KHL  model  were  detennined  for  all  three  alloys;  previously  detennined  values  were  used 
for  predictions  of  the  observed  material  response  in  strain-rate  jump  experiments  (all 
alloys)  and  high  strain  rate  experiments  at  different  temperatures  for  alloy  3  to  show  the 
applicability  of  the  model.  Two  of  the  six  constants  determined  were  kept  same  for  all 
three  alloys,  suggesting  that  the  strain  rate  sensitivity  and  the  strain  rate  effect  on  the 
work  hardening  of  the  alloys  were  independent  of  the  interstitial  content  and  small 
differences  in  the  initial  microstructures  of  these  alloys.  The  yield  strength,  work 
hardening  and  the  temperature  sensitivity  though  were  found  to  be  somewhat  different  for 
these  three  different  alloys  and  it  was  concluded  that  the  amount  of  interstitial  content 
and  the  initial  microstructure  had  been  responsible  for  these  variations  in  the  observed 
responses. 
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Table  1 


Chemical  composition  of  the  three  Ti-6A1-4V  alloys  used 


Material 

A1 

V 

Fe 

Y 

H 

N 

o 

C 

Ti 

Oeq 

Alloy  1 

6.26 

4.16 

0.14 

<0.0003 

0.0031 

0.008 

0.178 

0.047 

REM 

0.229 

Alloy  2 

6.30 

3.86 

0.18 

<0.0003 

0.0026 

0.003 

0.112 

0.045 

REM 

0.152 

Alloy  3 

5.97 

4.09 

0.15 

<0.0003 

0.0041 

0.008 

0.174 

0.043 

REM 

0.222 

Table  2 

KHL  model  material  constants  determined  for  the  three  Ti-6A1-4V  Alloys  used  in  the 
investigation 


A  B  ni  n0  C  in 

(MPa)  (MPa) 

1100  857.5  0.5455  0.6086  0.02204  1.6236 

988  747.1  0.5455  0.3986  0.02204  1.2214 

1069  874.8  0.5455  0.4987  0.02204  1.3916 


Ti-6A1-4V  Alloy  1 

Ti-6A1-4V  Alloy  2 
Ti-6A1-4V  Alloy  3 


Table  5 

Grain  Elongation  Ratio  and  average  grain  size  for  deformed  and  undeformed  samples  of 
alloy  3 


Sample 

Undeformed 

10'3  s'1  422  K 

Dynamic  296K 

Dynamic  422K 

Dynamic  715K 

Grain 

Elongation 

Ratio 

0.89 

0.36 

0.63 

0.50 

0.44 

Average 

Grain  Size 
(pm) 

9.86 

9.37 

9.80 

9.07 

7.91 

b) 


Figure  1.  (a)  Undeformed  and  (b)  18%  deformed  Ti-6A1-4V  alloy  2  at  296K  and  10°  s'1 
strain  rate  showing  the  globular  a  with  dispersed  (3  along  the  grain  boundaries. 


b) 


Figure  2.  (a)  Undeformed  and  (b)  20%  deformed  Ti-6A1-4V  alloy  3  at  422K  and  10"3  s'1 
strain  rate  showing  the  globular  a  with  dispersed  transformed  (3  (Widmanstatten 

structure). 


a) 


b) 


c) 


Figure  3.  Deformed  samples  of  alloy  3  at  high  strain  rates  and  different  temperatures,  a) 
at  296  K,  b)  at  422  K  and  c)  at  715  K  (compression  along  vertical  direction). 


Figure  4.  KFIL  model  correlations  with  the  room  temperature  (296K)  response  of  alloy  1 

at  different  strain  rates. 


Figure  5.  KHL  model  correlations  with  the  observed  response  of  alloy  1  at  different 
temperatures  at  constant  strain-rate  1CT4  s'1 
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Figure  6.  KHL  model  predictions  with  the  observed  response  of  alloy  1  during  a  strain- 
rate  jump  experiment  at  room  temperature  (296K).  Dynamic  predictions  are  adiabatic. 


Figure  7.  KHL  model  correlations  with  the  room  temperature  (296K)  response  for  alloy  2 

at  different  strain  rates. 


Figure  8.  KFIL  model  correlations  with  the  observed  response  of  alloy  2  at  different 
temperatures  at  constant  strain-rate  I O  4  s'1 


Figure  9.  KHL  model  predictions  with  the  observed  response  of  alloy  2  during  a  strain- 
rate  jump  experiment  at  room  temperature  (296K).  Dynamic  predictions  are  adiabatic. 
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Figure  10.  KHL  model  correlations  with  the  room  temperature  (296K)  response  for  alloy 
3  at  different  strain  rates.  (Published  data:  Khan  et  ah,  Int  J.  Plas  20,  2004) 


Figure  1 1 .  KHL  model  correlations  with  the  observed  response  of  alloy  3  at  different 

3  1 

temperatures  at  constant  strain-rate  10'  s'  .  (Published  data:  Khan  et  ah,  Int  J.  Plas  20, 
2004) 


Figure  12.  KHL  model  predictions  with  the  observed  response  of  alloy  3  during  a  strain- 
rate  jump  experiment  at  room  temperature  (296K).  Dynamic  predictions  are  adiabatic. 
(Published  data:  Khan  et  ah,  Int  J.  Plas  20,  2004). 


Figure  13.  KHL  model  predictions  with  the  observed  response  of  alloy  3  at  high  strain 
rates  and  different  temperatures.  All  the  predictions  are  adiabatic. 
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Abstract 

Multiaxial  experiments  on  Ti-6A1-4V  titanium  alloys  are  presented.  Different  loading 
conditions  are  applied  in  order  to  get  the  comprehensive  response  of  the  alloy.  The  strain 
rates  are  varied  from  the  quasi-static  regimes  to  dynamic  regimes  and  the  corresponding 
material  response  is  obtained.  The  specimen  is  deformed  to  large  strains  in  order  to  study 
the  material  behavior  under  finite  deformation  at  various  strain  rates.  Torsional  Kolsky 
bar  is  used  to  achieve  material  strain  rates  up  to  1000  s'1.  The  non-proportional  loading 
experiments  comprise  of  an  initial  uniaxial  loading  to  a  certain  level  of  strain  then 
followed  by  biaxial  loading,  using  a  channel-type  die  at  various  rates  of  loadings,  as  well 
as  dynamic  torsion  followed  by  dynamic  compression.  All  the  non-proportional 
experiments  are  carried  out  at  room  temperature.  Experiments  are  also  performed  to 
investigate  the  anisotropic  behavior  of  the  alloy.  An  orthotropic  yield  criterion  proposed 
by  Cazacu  et  al.  (2005)  for  anisotropic  hexagonal  closed  packed  materials  with  strength 
differential  is  used  to  generate  the  yield  surface.  Based  on  the  definition  of  the  effective 
stress  of  this  yield  criterion,  the  observed  material  response  for  the  different  loading 
conditions  under  large  deformation  is  modeled  using  the  Khan-FIuang-Liang  (KHL) 
equation  assuming  isotropic  hardening.  The  model  constants  used  were  pre-determined 
from  the  extensive  uniaxial  experiments  presented  in  the  earlier  paper  [Int.  J.  Plast.,  20 
(2004),  2233-2248].  The  model  predictions  are  found  to  be  extremely  close  to  the 
observed  material  response. 
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1.  Introduction 


Titanium  alloys  have  found  a  wide  variety  of  usage  in  the  aerospace  industry  because 
of  its  high  strength  to  weight  ratio.  Medical  industry  has  uses  for  titanium  in  surgical 
equipment  and  implants  because  of  its  biocompatibility.  Petrochemical  industry  has  also 
found  immense  usage  for  this  alloy  because  of  its  high  corrosion  resistance.  The  defense 
industry  is  interested  in  advancing  the  application  of  these  alloys  to  armor  systems.  This 
is  because  of  the  high  strength  to  weight  ratio  and  the  corrosion  resistance  of  these  alloys 
encourage  its  use  in  armor  of  tanks  and  as  compartments  for  the  fuel  and  ammunitions. 
These  alloys  come  in  different  grades  depending  on  the  impurity  contents.  The  ELIs  or 
the  extra  low  interstitials  are  more  popular  in  ground  applications  because  of  their 
increased  deformability  and  fracture  resistance  thus  favoring  the  ballistic  properties. 
Although  this  grade  of  alloy  has  a  slightly  lesser  strength  than  the  ones  used  for 
aerospace  applications,  however,  it  is  cheaper  than  conventional  Ti-6A1-4V  titanium 
alloy. 


A  large  number  of  studies  have  been  carried  out  in  understanding  and  modeling  the 
response  of  these  titanium  alloys  (Follansbee  and  Gray,  1989;  Lesuer,  2000;  Nemat- 
Nasser  et  al.  2001;  Picu  and  Majorell,  2002;  Lee  and  Lin,  1997;  Khan  et  al.  2004).  These 
studies  were  focused  mainly  on  the  uniaxial  response  of  the  alloys  at  different  strain  rates 
and  temperatures.  However,  it  is  well  known  that  the  stress  states  in  these  alloys  while  in 
use  are  definitely  more  complex  and  multiaxial  in  nature.  Hence  it  becomes  extremely 
important  to  understand  the  material  behavior  during  these  loading  conditions.  The 
uniaxial  material  responses  and  modeling  act  only  as  the  first  step  towards  understanding 
the  actual  more  complex  conditions  that  exist  and  how  they  affect  the  corresponding 
responses. 

There  have  been  a  few  studies  understanding  the  torsional  behavior  of  the  Ti-6A1- 
4V  titanium  alloys  which  are  available  in  the  literature.  Liao  and  Duffy  (1998), 
performed  a  series  of  high  strain  rate  (y~10  /sec)  torsion  experiments  to  study  the 
process  of  initiation  and  formation  of  adiabatic  shear  bands  in  conventional  Ti-6A1-4V 
titanium  alloy.  The  temperature  rise  in  the  alloys  during  deformation  was  also  measured. 
Macdougall  and  Harding  (1999)  have  also  performed  a  series  of  torsional  experiments  on 
thin  tubular  Ti-6A1-4V  specimens  to  achieve  the  response  of  these  alloys  in  shear.  They 
also  measured  the  temperature  rise  in  the  specimens.  Constitutive  modeling  was 
performed  using  Zerilli-Armstrong  (ZA)  model;  high  strain  rate  torsional  response  of  the 
alloys  and  the  material  constants  were  incorporated  into  the  ABAQUS /explicit  FE  code  to 
predict  the  tensile  response  of  the  alloys.  They  concluded  that  the  ZA  model  was  unable 
to  model  the  markedly  reduced  work  hardening  behavior  at  high  strain  rate  torsion 
deformation.  They  also  suggested  that  strain  rate  and  temperature  dependence  of  work 
hardening  required  a  modification  in  the  model.  Chichili  et  al.  (2004)  have  studied  the 
torsional  high  strain  rate  response  of  a-titanium.  They  perfonned  compression-torsion 
Kolsky  bar  experiments  at  high  rates.  The  specimens  used  in  their  study  were  of  two 
different  kinds.  First  one  was  the  conventional  thin  tubular  specimen  (Marchand  and 
Duffy,  1988)  and  the  other  specimen  is  solid  and  had  a  circumferential  notch  in  the  test 
section  to  study  shear  localization. 


Non-proportional  loading  experiments  are  useful  in  determining  the  multiaxial 
response  of  the  material  when  subject  to  complex  loading  conditions.  Bridgman  (1946) 
was  first  to  study  the  material  response  during  two  dimensional  compressive  experiments 
on  steel.  He  used  two  independent  hydraulic  presses  to  apply  compressive  loads  to  two 
mutually  perpendicular  faces  of  a  rectangular  block  specimen,  the  third  face  was  free  to 
deform.  Florenz  (1969)  modified  this  experiment  to  replace  one  of  the  hydraulic  presses 
with  almost  rigid  tool  steel  walls.  The  load  was  applied  in  one  direction  and  the  tool  steel 
walls  applied  a  stress  on  that  face  by  constraining  the  deformation  in  that  direction; 
however  he  did  not  measure  the  stress  in  the  constrained  direction.  This  version  of 
experiments  is  often  called  the  channel  die  experiments.  This  experiment  was  modified 
by  Khan  and  Wang  (1990)  in  which  all  the  non-zero  stresses  and  strains  were  accurately 
measured  in  all  three  directions  for  annealed  aluminum.  Later  it  was  used  by  Khan  and 
Liang,  (2000)  to  observe  the  response  of  Aermet,  Ta  and  Ta-2.5W  under  these  loadings. 

Titanium  alloys  in  the  low  and  medium  temperature  regime  consist  mainly  of  the 
HCP  a  -phase  with  very  little  dispersed  (BCC)  /?  -phase  in  between  the  equiaxed  a 
grains.  These  hexagonal  closed  packed  (HCP)  metals  are  known  to  display  plastic 
anisotropy  and  a  strong  strength  differential  in  tension  and  compression.  This  is  because 
of  the  interaction  between  crystallographic  slip  and  deformation  twinning  (Cazacu  et  ah, 
2005).  Twinning  is  also  known  to  be  dependent  on  the  sign  of  the  shear  stresses 
(Hosford,  1998).  Also,  anisotropy  exists  in  these  alloys  mainly  due  to  the  strong 
crystallographic  texture  induced  due  to  the  rolling  process  (Hosford,  2005).  Therefore, 
multiaxial  modeling  of  these  hep  titanium  alloys  cannot  be  performed  using  definitions  of 
the  effective  stress,  strain  and  strain  rate  based  on  the  isotropic  Jj-flow  theory  as  these 
alloys  have  been  known  to  be  anisotropic  in  nature.  Similar  conclusions  were  found  by 
Chichili  et  al.  (2004)  while  studying  the  multiaxial  response  of  a  -titanium  . 

Anisotropy  has  been  studied  in  great  detail  by  many  researchers  over  the  past  many 
decades.  Hill  (1948)  was  the  first  to  quantitatively  determine  the  yield  surface  for 
anisotropic  materials.  Hill’s  anisotropic  model  was  primarily  based  on  the  generalization 
of  the  von-Mises  yield  criterion.  Although  this  model  has  been  useful  because  of  its 
simplicity  it  was  not  able  to  capture  the  observed  anisotropy  in  yield  stress  and  the  R- 
ratios  together  (Lademo  et  al.  1999).  Other  prominent  anisotropic  yield  criteria  are  by 
Hosford  (1979),  Barlat  and  co  workers  (Barlat  et  al.  1991,  1997,  2003),  Karafrllis  and 
Boyce  (1993),  Yoon  et  al.  (2006),  etc.  These  criteria  dealt  with  the  modeling  of  cubic 
crystals  and  assume  same  yield  stress  in  tension  an  compression.  Some  anisotropic 
models  have  been  proposed  and  used  in  the  literature  for  hep  polycrystals  by  Tome  and 
Lebensohn  (2004)  and  Staroselsky  and  Anand  (2003).  Macroscopic  criteria  are  generally 
less  time  consuming,  such  as  the  orthotropic  criterion  by  Cazacu  et  al.  (2005)  which  is 
based  on  the  linear  transformation  of  the  Cauchy  stress  tensor  deviator. 

The  present  study  includes  a  comprehensive  set  of  experiments  for  understanding 
the  response  of  this  alloy  and  consequently  modeling  the  observed  responses  of  these 
alloys  at  different  strain  rates.  It  is  an  extension  of  the  previous  investigation  by  Khan  et 
al.  (2004).  In  this  investigation,  the  high  strain  rate  torsion  experiments  were  carried  out 
using  the  modified  Kolsky  bar  apparatus  (Hartley  et  al.,  1985).  This  was  followed  by 


compression  on  the  same  samples  again  at  high  strain  rates  albeit  at  different  operating 
temperatures.  The  quasi-static  torsion  experiments  were  carried  out  using  the  MTS  axial/ 
torsion  material  testing  system.  The  non-proportional  loading  experiments  were 
performed  using  the  channel  die  at  small  and  intermediate  strain-rate  regime.  During 
these  experiments  the  stress  state  was  designed  to  change  from  uniaxial  to  biaxial 
compression  and  the  corresponding  response  was  observed  and  recorded.  Constitutive 
modeling  of  these  observed  responses  was  performed  using  the  KHL  model.  The 
experimental  data  was  first  converted  into  effective  stress  and  effective  strain  based  on 
definitions  of  the  orthotropic  yield  criteria  proposed  by  Cazacu  et  al  (2005).  For  an 
alternative  constitutive  modeling  technique  the  reader  is  referred  to  Abed  and  Voyiadjis, 
(2005). 

2.  Experimental  Procedure 

2. 1  Material 

The  titanium  alloy  used  in  this  study  is  Ti-6A1-4V  with  high  oxygen  content.  The 
equivalent  oxygen  content,  O  =  O+2N+0.75C  (DeMeester  et  al.  1975),  of  this  alloy  is 

0.222  %.  This  alloy  was  manufactured  using  single  electron  beam  melting  process, 
resulting  in  higher  oxygen  content  than  the  ELI  grade.  The  chemical  composition  of  this 
titanium  alloy  is  shown  in  Table  1. 

2.2  Quasi-static  free-end  torsional  experiments 

Torsion  specimens  for  quasi-static  loading  were  machined  as  shown  in  Figure  1(a). 
The  specimens  were  in  the  form  of  a  thin  hollow  tube  with  a  square  flange  to  minimize 
any  radial  variation  in  the  shear  stress  in  the  test  section  (Hosford,  2005).  A  female  slot 
similar  to  the  square  flange  was  machined  in  the  fixture  holding  the  specimen.  The 
geometry  used  is  somewhat  similar  to  the  hexagonal  flanges  used  by  Liao  and  Duffy, 
(1998)  for  dynamic  torsional  experiments.  Strain  gage  rosettes  were  mounted  on  the  test 
section  of  the  specimen.  These  gages  were  bonded  at  an  angle  of  ±45°  degrees  to  the 
vertical  axis.  The  torque  experienced  by  the  torsional  transducer  during  rotation  of  the 
lower  fixture  in  an  angle  controlled  mode,  and  the  axial  displacement  were  also  recorded. 
The  shear  stress  and  strain  were  calculated  as  in  Khan  and  Liang,  (2000).  The  axial  force 
was  zero  during  the  experiment.  Experiments  were  performed  at  different  strain  rate 
using  this  experimental  setup.  The  gage  section  wall  thickness  and  length  were  0.018  in. 
(0.4752  mm)  and  0.125  in.  (3.175  mm),  respectively.  In  one  experiment  the  shear  strain 
rate  was  kept  constant  at  1.732xl0~5  s'1.  A  strain  rate  jump  experiment  was  also 
performed  where  the  shear  strain  rate  was  abruptly  changed  from  initially  1.732 xlO’V1 
to  1.732  s'1  at  approximately  8.5  %  plastic  shear  strain. 

2.3  Dynamic  torsional  experiments 

Dynamic  torsion  experiments  were  performed  using  the  modified  torsional  Kolsky 
bar  (TSHB)  setup.  The  specimens  were  machined  as  shown  in  Figure  1(b).  The  schematic 
diagram  of  the  apparatus  used  is  also  shown  in  Figure  2.  The  torsional  Kolsky  bar  works 


on  the  same  principle  as  the  compression  bar  with  the  following  differences.  A  section  of 
the  incident  bar  CD  (Figure  2)  is  quasi-statically  pre-torqued.  The  bolt  holding  the  clamp 
at  D  is  ruptured,  thus  releasing  a  shear  wave  in  the  incident  bar.  The  strain  gages  on  the 
incident  and  transmitted  bars  are  mounted  at  ±  45°  to  the  axis  of  these  bars.  These  strain 
gages  are  connected  to  a  Nicolet  440  digital  oscilloscope  using  a  potentiometer  circuit. 
The  specimen  is  glued  to  the  incident  and  transmitted  bars  using  a  high  strength  adhesive. 
The  shear  stress  and  strain  rate  are  obtained  by  using  Eqns.  (1)  and  (2). 

T{t)=G^i:yT  -(1);  r(0  =  7^r«  --(2) 

d;h  LRh 

where  G  is  the  shear  modulus  of  the  bar,  f?/,  is  the  bar  diameter,  c'0  is  the  shear  wave 
speed,  ds  is  the  mean  diameter  of  the  specimen,  L  and  h  are  the  length  and  the  thickness 
of  the  test  section  of  the  specimen,  respectively.  Also,  Yt  and  Yr  are  the  transmitted  and 
reflected  shear  strains,  respectively.  The  shear  strain  rate  is  integrated  over  time  to  get  the 
shear  strains.  Correlating  the  stress  and  strain  with  respect  to  time  gives  the  shear  stress- 
strain  response  of  the  material  under  dynamic  torsional  loading. 

The  diameter  of  the  incident  and  transmitted  bars  used  in  the  experiment  is  0.84  in 
(21.34  mm).  These  bars  were  made  from  hardened  Vascomax-C350  steel.  Hollow  tubular 
specimens  are  used  as  in  the  case  of  quasi-static  torsional  experiments.  The  inside 
diameter  of  the  specimen  is  approximately  0.5  in.  (12.7  mm)  and  the  outside  diameter  is 
0.536  in.  (13.6  mm),  approximately.  A  more  detailed  discussion  of  this  experimental 
technique  can  be  found  in  the  literature  (Hartley  et  al.  1985). 

2.4  Dynamic  torsion-compression  experiments 

Dynamic  torsion  compression  experiments  were  performed  on  the  titanium  alloys 
using  the  TSHB  technique  as  described  above.  Once  the  hollow  tubular  specimen  was 
deformed  plastically  under  torsional  loading,  the  test  section  was  machined  in  a  form  of  a 
circular  ring.  The  ring  was  used  as  a  specimen  and  compressed  under  dynamic  loading 
conditions  using  the  compression  split-Hopkinson  bar  technique  (SPHB).  This  technique 
was  used  for  obtaining  the  response  of  the  specimens  at  room  temperature  and  also  at 
different  temperatures.  The  details  of  the  SPHB  technique  can  be  found  in  Liang  and 
Khan  (1999)  and  also  in  Follansbee  (1979).  All  dynamic  torsion  experiments  were 
performed  at  room  temperature;  however,  the  subsequent  dynamic  compression 
experiments  using  the  rings  machined  from  the  specimens  after  dynamic  torsional 
loadings  were  performed  at  room  and  higher  temperatures.  In  these  experiments,  the 
specimens  were  heated  to  known  temperatures  (3  OOF  and  600F)  and  then  loaded  in 
compression  at  dynamic  strain  rates.  The  diameter  of  the  incident,  transmitted  and 
projectile  bars  of  the  split-Hopkinson  pressure  bar  (SHPB)  apparatus  was  0.5  in.  (12.7 
mm).  These  bars  were  made  from  Vascomax-C350  steel  heat  treated  to  maximum 
hardness. 


2.5  Non-Proportional  Loading  Experiments 


Non-proportional  compressive  loading  experiments  were  performed  on  the  titanium 
alloys  at  strain  rates  of  10'  ,  10"  and  1  s'  in  the  loading  direction.  The  specimens  were 
machined  in  the  form  of  a  rectangular  block  of  dimensions  0.72  in.  (18.29  mm)  in  height 
(loading  direction),  0.4  in.  (10.16  mm)  and  0.5  in.  (12.7  mm)  in  the  other  two  directions. 
The  schematic  diagrams  of  the  channel  die  for  the  non-proportional  loading  experiments 
and  the  rectangular  block  specimen  are  shown  in  Figure  3.  It  consists  of  a  hardened 
Vascomax-C350  steel  die  with  removable  side  plate  inserts,  a  removable  base  plate  and  a 
plunger  or  loading  block  also  made  from  the  same  material.  The  plate  inserts  and  base 
plate  are  made  removable  so  that  the  specimen  can  be  taken  out  easily  after  the 
experiment.  The  specimen  is  placed  inside  the  die  as  shown  in  the  apparatus.  Initially  the 
specimen  is  subjected  to  uniaxial  deformation  in  direction  1  (before  the  gap  between 
specimen  and  die  plate  inserts  in  direction  2  closes).  The  direction  1  is  the  loading 
direction  as  shown  in  the  schematic.  Then,  after  some  pre-designed  finite  deformation, 
the  block  is  subjected  to  large  deformations  under  biaxial  compression  due  to  constraint 
in  the  other  direction  (direction  2).  The  3rd  direction  is  free  of  stress  as  the  specimen 
deforms  freely  in  that  direction.  The  amount  of  finite  uniaxial  deformation  can  be 
controlled  by  choosing  the  right  dimension  of  the  specimen  in  the  2nd  direction 
(approximately  10%  in  the  present  experiments).  Uniaxial  loading  is  applied  through  the 
MTS  material  testing  system  through  a  hardened  loading  block  also  called  the  plunger. 
The  dimension  of  this  plunger  is  selected  so  that  there  is  neither  any  pinching  nor  any 
appreciable  plastic  flow  through  any  space  between  the  plunger  and  the  hardened  steel 
plates.  Strain  gages  were  mounted  on  the  specimen  in  the  loading  direction  to  record  the 
strain  in  that  direction.  Clip  gages  were  used  to  record  the  strain  in  the  direction  3,  or  the 
direction  in  which  the  specimen  is  free  to  defonn.  These  gages  are  made  by  bonding 
uniaxial  strain  gages  at  exactly  the  same  location  on  the  tension  and  compression  side  of 
a  cantilever  beam  made  from  very  thin  blade  strips  of  high  strength  steel.  These  strips  are 
connected  to  the  die  in  the  form  of  a  cantilever  beam  as  shown  in  Figure  3.  Steel  pins 
were  used  to  connect  the  specimen  to  the  blade  strips  of  the  clip  gages.  Once  the  load  is 
applied  on  the  specimen  in  the  direction  1 ,  the  pins  start  to  move  outward  in  the  direction 
3  due  to  the  Poisson  effect.  This  causes  the  clip  gages  or  the  cantilever  beams  (blade 
strips)  to  bend  outward.  The  gages  on  the  inside  surface  of  the  beams  (or  strips) 
experience  tensile  elastic  deformation  while  the  gages  on  the  outside  surface  experiences 
an  equal  amount  of  compressive  strains.  The  gages  are  then  connected  in  a  full 
Wheatstone  bridge  arrangement  so  the  magnitude  of  these  strains  are  added  and  thus 
quadrupled.  These  cantilever  beam  clip  gages  are  calibrated  before  the  experiment  by 
moving  them  (known  amounts)  to  get  strains  in  the  specimen  in  this  direction.  In  order  to 
measure  the  stress  in  the  direction  2,  high  pressure  manganin  gages  (LM-SS-1 10FB-048) 
manufactured  by  Vishay  Micro-Measurements  Division  are  used.  These  gages  are 
bonded  to  the  surface  of  the  removable  steel  plate  inserts,  facing  the  die.  Once  the  level 
of  deformation  in  the  specimen  is  such  that  it  starts  to  push  against  the  hardened  steel 
plates  these  gages  experience  an  increase  in  the  pressure  (stress),  which  is  recorded  using 
the  teststar  software  of  the  MTS.  These  gages  are  also  calibrated  prior  to  experiment.  The 
amount  of  lubrication  is  very  critical  on  the  removable  walls,  the  base  plate  and  plunger 
during  these  experiments,  as  experienced  by  others  as  well  (Bridgman,  1946  and  Bell, 


1988).  In  our  case,  high  vacuum  grease  manufactured  by  Dow  Corning  was  used  together 
with  graphite  powder.  This  ensured  that  there  was  no  bulging  or  reverse  barreling  during 
the  experiment  and  thus  the  deformation  was  uniform. 

2.6  Uniaxial  Experiments  in  Different  Directions 

In  order  to  determine  the  level  of  anisotropy  in  the  alloy,  uniaxial  quasi-static 
compression  experiments  were  performed  in  different  directions  of  the  plate.  The 
principal  axes  of  orthotropy  are  denoted  x,  y  and  z,  where  x  axis  direction  represents  the 
rolling  direction  of  the  plate,  y  direction  is  transverse  to  the  rolling  direction  in  plane  of 
rolling  and  z  direction  is  the  thickness  direction  of  the  plate.  For  the  determination  of  the 
in-plane  anisotropy,  uniaxial  quasi-static  compression  experiments  were  performed  in  the 
rolling  direction,  45°  to  rolling  direction  and  transverse  to  rolling  direction.  Strain  values 
were  recorded  for  the  loading  and  transverse  to  loading  directions  also  in  order  to 
calculate  the  corresponding  R-ratios.  Additionally,  tension  experiments  were  also 
performed  in  the  rolling  and  transverse  to  rolling  directions.  All  the  experiments  for 
anisotropic  analysis  were  performed  at  a  constant  strain  rate  of  10'  s'  .  The  0.2%  offset 
definition  of  yield  was  used  for  calculation  of  the  yield  strength  of  the  alloy.  This 
definition  is  found  suitable  for  this  alloy,  as  unlike  FCC  metals  (aluminum,  etc)  it  has  a 
well  defined  yield  stress  at  0.2  %  offset. 

3.  Orthotropic  Yield  Criteria  and  Constants 

In  order  to  describe  the  strength  differential  and  the  anisotropy  typically  shown  by 
hexagonally  closed  packed  (hep)  metals  and  alloys,  Cazacu  et  al.  (2005),  proposed  a  new 
phenomenological  orthotropic  criterion  that  can  capture  the  response  of  the  material  in 
different  directions  accurately.  The  authors  developed  the  criterion  initially  for  materials 
with  strength  differential,  i.e.,  different  strength  in  tension  and  compression  and  extended 
it  for  anisotropic  materials.  This  criterion  is  based  on  the  linear  transformation  of  the 
deviatoric  part  of  the  Cauchy  stress  tensor,  S  similar  to  previous  studies  by  Barlat  and  co¬ 
workers  (Barlat  et  al.  199 1 ;  Barlat  et  al.  1997;  etc)  and  has  previously  been  used  by  many 
researchers  (Lademo  et  al.,  1999;  Lademo  et  al.,  2002;  Abedrabbo  et  al.,  2006,  etc) 
successfully.  The  transformed  tensor  X  is  defined  as: 

I  =  C[S]  (3) 

where  C  is  a  4  th  order  tensor.  The  ortho  tropic  yield  criterion  is  given  as: 

(Z,  I  ■ -  k.  Z, )"  +  Qz,  I  -  k.  Z2 )“  +  Qz,  |  -  k.  I , )"  =  F  (4) 

The  exponent  a,  in  the  yield  criterion  is  considered  to  be  a  positive  integer  (in  our 
case  a= 2).  The  parameter  k  is  a  material  constant,  which  for  a  fixed  a,  is  determined  by 
the  ratio  of  the  strength  in  tension  and  compression.  The  function  F,  gives  the  size  of  the 
yield  locus.  F  can  be  defined  either  as  a  constant  or  a  function  of  the  total  plastic  strain 
(isotropic  hardening).  The  parameter  k,  is  defined  as; 

1-/7 


where  h  is  given  by, 


For  a  2-D  case,  the  transformation  tensor  C,  is  given  by 
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The  principal  anisotropic  stress  tensor  components  are  given  as: 


(10) 


According  to  the  linear  transformation,!]  =  C[s],  the  non-zero  components  of  the 
anisotropic  stress  tensor  are  given  as, 


Where,  sx,  sy,  sxy  and  sz  are  the  components  of  the  deviatoric  Cauchy  stress  tensor.  Upon 
expansion  Eqn  (11)  leads  to 
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Applying  the  various  stress  states,  the  yield  stress  in  tension  and  compression  along  the 
rolling  direction  x  (0°),  according  to  the  orthotropic  criterion  are  given  as 
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where, 
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Similarly,  the  tensile  and  compressive  yield  stresses  in  the  transverse  direction  (90°),  y, 
are  given  by, 
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The  stress  state  in  material  coordinates  for  a  compression  experiment  in  the  45°  direction 

c  c  c 

V  45  ^45  i  ^  45 

is  given  as  =  -  — ,  =  -  —  and  crvv  =  -  — 

Based  on  this,  the  principal  values  of  the  transformed  tensor  for  this  experiment  are  given 
as: 
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where, 
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Finally,  the  yield  stress  for  compression  experiment  at  45°  to  rolling  direction  is  given  as 
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Furthermore,  yielding  under  pure  shear  parallel  to  the  orthotropy  axes  is  given  by 
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The  R-ratio  for  different  orientations  of  the  specimen  is  defined  as  the  ratio  of  the 
incremental  transverse  to  thickness  strain  and  is  given  as 

(27) 
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de. 


Through  constancy  of  volume  assumption  this  can  also  written  as 
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here,  dewandde,  are  the  incremental  width  and  longitudinal  plastic  strain  respectively. 

For  isotropic  materials  the  R-ratios  are  normally  equal  to  unity.  Through  intense 
mathematics  and  use  of  associated  flow  rule  (incorporating  the  assumption  that  the  plastic 
potential  coincides  with  the  yield  function),  the  R-ratios  can  be  written  for  the  different 
directions  as  follows 
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The  exponent  a,  is  chosen  to  be  2,  as  it  was  also  used  by  Cazacu  et  al.  (2005)  for 
other  Ti  alloys.  This  orthotropic  yield  criterion,  for  a  fixed  exponent  a,  utilizes  7  material 
constants  for  a  plane  stress  (2D)  case.  Correspondingly,  detennination  of  these  material 
constants  requires  seven  experimental  results.  For  a  more  detailed  explanation  of  this 
yield  criterion,  the  reader  is  referred  to  the  original  paper  by  Cazacu  et  al.  (2005). 

Experimental  results  from  uniaxial  tension  and  compression  yield  stresses  in  rolling 
and  transverse  to  rolling  directions,  uniaxial  compression  in  the  45°  to  the  rolling 
direction  and  R-ratios  in  the  rolling,  and  transverse  to  rolling  directions  for  the 
compression  case  are  used.  These  material  constants  were  calculated  using  least  square 
optimization  technique  in  MATLAB.  The  value  of  the  material  constant  k,  used  in  this 
investigation  was  calculated  from  the  ratio  of  the  yield  in  tension  and  compression  in  the 
transverse  to  rolling  direction.  Table  2,  shows  the  material  constants  which  were 
calculated  using  least  square  optimization  method.  The  yield  stresses  determined  from 
these  experiments  and  R-ratios  along  with  the  values  predicted  by  the  model  based  on  the 
optimized  material  constants  can  be  seen  in  Table  3.  Figure  4  shows  the  yield  surface  for 
the  orthotropic  yield  criteria  based  on  the  material  constants  calculated  for  our  case. 
Figure  5  shows  the  normalized  yield  stress  correlations  from  the  orthotropic  criterion 
along  with  the  normalized  experimental  yield  stresses  at  in  different  directions  in  plane  of 
the  plate.  It  can  be  seen  from  the  Figure  4  and  5  and  Table  3  that  the  orthotropic  model 
works  reasonably  well  for  hep  metals  in  the  development  of  the  yield  surface  and  that  the 
values  of  R-ratios  (which  incorporates  the  flow  rule)  correlated  by  this  model  are 
extremely  close  to  the  experimental  data.  This  implies  that  it  is  reasonable  to  use 
effective  stress,  strain  and  strain  rate  definitions  using  this  ortho  tropic  model. 


4.  Experimental  Results  and  Constitutive  Model  Predictions 


The  Khan-Huang-Liang  constitutive  model  is  used  to  simulate  the  observed 
behavior  of  the  titanium  alloys  under  the  complex  multiaxial  loadings  at  various  strain 
rates  and  temperatures.  The  detailed  modeling  procedure  for  the  above  model  can  be 
found  in  Khan  and  Liang,  (1999)  &  Khan  et  al.  (2004).  The  constitutive  model  is 
expressed  as 
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where,  <J,£P,  £  are  the  effective  stress,  effective  plastic  strain  and  strain  rate, 
respectively.  Tm,  T,  Tref  are  melting,  current,  and  reference  temperatures  in  degrees 

Kelvin,  respectively.  D%= 106  s' ’(arbitrarily  chosen  upper  bound  strain  rate)  andf*=l  s'1. 
A,B,nvn0,C  and  m  are  the  material  constants.  For  Ti-6A1-4V  alloys,  the  melting 

temperature  was  taken  to  be  1933  K  (ASM  Handbook,  1994).  The  reference  temperature 
was  room  temperature  (296K).  These  material  constants  have  been  previously  determined 
from  uniaxial  experimental  results  at  different  strain  rates  and  temperatures  for  the  same 
alloy  (Khan  et  al.,  2004). 


The  experimental  results  from  quasi-static  mono  tonic  torsion,  quasi-static  jump 
torsion  and  dynamic  torsion  experiments  are  shown  in  Figure  6,  along  with  the  model 
predictions  (the  material  constants  were  not  changed  from  the  published  values,  Khan  et 
ah,  2004).  The  experimental  results  in  these  multiaxial  stress  states  were  converted  into 
effective  stress  ( a e )  using  the  definition  of  the  yield  function  (Eqn.  4),  so  that  the 
complexity  of  the  stress  state  reduces  to  a  single  scalar  quantity. 

In  a  similar  manner  the  effective  strain,  ( £e ),  and  strain  rates,  (£e),  are  also 

calculated  from  the  experimental  results  by  considering  the  principle  of  incremental  work 
per  unit  volume  (Khan  and  Huang,  1995).  The  experimental  results  were  then  plotted 
against  the  model  predictions  based  on  these  definitions. 

dWp  =  (Jydefj  =  crjd£p  =  aed£p  (28) 

The  KHL  model  is  able  to  predict  the  response  of  these  quasi  static  as  well  as  the 
dynamic  torsion  experiments  reasonably  well  (Fig.  6).  The  KHL  predictions  are  made 
assuming  adiabatic  defonnations  for  the  1.732  s-1  (shear  strain  rate)  as  well  as  the 
dynamic  experimental  results  as  an  increase  in  the  temperature  were  noticed  during  both 
experiments.  The  reader  is  referred  to  another  paper  (Khan  et  ah,  2004)  for  a  detailed 
discussion  of  the  calculation  of  the  temperature  rise  in  the  specimen. 

Figure  7  shows  the  observed  response  of  the  alloy  under  dynamic  torsion  loading 
followed  by  dynamic  compression  at  room  temperature.  Here  again,  the  stress  obtained 
using  the  KHL  constitutive  model  is  converted  for  adiabatic  case  for  both  the  loading 
conditions  (torsion  and  compression)  as  the  effective  strain  rate  of  the  experiment  is  in 
the  dynamic  regime  and  thus  the  deformation  is  adiabatic  resulting  in  thermal  softening. 
The  KHL  model  predictions  are  again  found  to  be  in  excellent  agreement  with 
experimentally  observed  response  of  the  alloy. 

Figure  8  shows  the  measured  response  of  the  alloy  under  dynamic  torsion  followed 
by  dynamic  compression  at  higher  temperature  (422  K).  The  KHL  model  predictions 
were  converted  for  adiabatic  case  as  is  done  in  the  previous  case.  The  predictions  are 
again  in  very  good  agreement  with  the  actual  observed  material  response  for  both  types 
of  dynamic  loading  conditions  showing  the  capability  of  the  model  in  capturing  the 
observed  behavior  of  the  alloy  under  complex  loading  conditions  when  used  in 
conjunction  with  the  anisotropic  criterion. 

Figure  9  shows  the  experimental  results  of  the  alloy  during  dynamic  torsion 
followed  by  dynamic  compression  at  588  K.  The  model  predictions  are  again  converted 
to  adiabatic  case  for  both  loadings.  The  differences  between  predictions  and  the  observed 
responses  are  reasonably  small. 

The  true  stress  and  strain  response  of  the  alloy  under  non-proportional  loading 
conditions  at  I  0"1  s'1  strain  rate  in  the  loading  direction  is  shown  in  Figure  10.  The  rise  in 
the  stress  in  the  loading  direction  is  clearly  visible  as  the  specimen  becomes  constrained 
once  a  certain  predetermined  level  of  deformation  in  the  loading  direction  is  achieved. 


Consequently,  there  is  a  corresponding  rise  in  the  stress  in  the  constrained  direction  until 
the  end  of  the  experiment.  Figure  1 1  shows  the  effective  stress  and  effective  strain 
responses  for  different  strain  rates  in  the  loading  directions  along  with  the  model 
predictions.  Effective  stresses,  strains  and  strain  rates  are  calculated  based  on  the 
orthotropic  yield  criterion  of  Cazacu  et  al  (2005)  and  the  principle  of  incremental  work 
per  unit  volume.  The  stresses  and  strains  are  then  plotted  with  the  model  predictions. 
Again,  the  published  values  of  constants  [Khan  et  al  (2004)]  are  used  and  are  not 
modified  using  the  present  experimental  results.  The  model  prediction  for  the  case  of  10'1 
s'1  strain  rate  in  the  loading  direction  is  converted  assuming  adiabatic  deformation,  as 
there  was  a  significant  increase  in  the  temperature  in  the  specimen  noticed  during  the 
deformation.  The  model  predictions  are  in  excellent  agreement  with  the  experimentally 
observed  responses  for  all  strain  rates.  The  strain  rates  used  in  the  model  calculations  in 
all  non  proportional  loading  experiments  are  effective  strain  rates,  and  are  computed 
using  the  three  strain  histories  along  with  the  principle  of  incremental  work  per  unit 
volume. 

5.  Conclusions 

Quasi  static  and  dynamic  torsion  experimental  results,  dynamic  torsion  loadings 
followed  by  dynamic  compression  experiments  at  different  strain  rates  and  temperatures 
were  presented  in  this  study  on  a  single  electron  beam  melt  Ti-6A1-4V  titanium  alloy. 
Non  proportional  loading  experiments  at  various  strain  rates  were  also  presented.  These 
experimental  results  formed  a  comprehensive  data  to  display  the  behavior  of  the  alloy 
under  various  complex  loading  conditions. 

Using  the  set  of  model  constants  for  this  alloy  which  were  calculated  from  the 
uniaxial  responses  at  different  strain  rates  and  temperatures  in  an  earlier  paper  (Khan  et 
al.  2004),  the  capability  of  the  model  to  predict  accurately  the  observed  new  responses 
under  the  complex  multiaxial  loading  conditions  was  demonstrated.  The  model  was  able 
to  predict  the  material  response  under  pure  torsion  and  torsion  followed  by  compression 
extremely  well  using  the  effective  stress,  strain  and  strain  rate  definitions  from  the 
ortho  tropic  yield  criteria  and  principle  of  incremental  work  per  unit  volume.  In  the  non 
proportional  channel  die  experiments,  the  model  was  again  able  to  predict  the  observed 
material  responses  under  non-proportional  loading  conditions  reasonably  well. 
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Table  1 


Chemical  composition  of  the  Ti-6A1-4V  alloy  used 


Material 

A1 

V 

Fe 

Y 

H 

N 

o 

C 

Ti 

Ti-6A1-4V 

Titanium  alloy 

5.97 

4.09 

0.15 

<0.0003 

0.0041 

0.008 

0.174 

0.043 

REM 

Table  2 


Material  constants  for  the  titanium  alloy  using  the  orthotropic  yield  criteria  of  Cazacu  et 
ah, (2005) _ 


a 

k 

Cn 

C12 

C13 

C22 

c23 

c33 

c66 

2 

-0.0513 

1.5717 

2.7749 

1.118 

1.6907 

1.2405 

0.5651 

1.2133 

Table  3 

Experimental  and  Predicted  values  of  the  yield  stresses  and  R-ratios  under  different 


loading  conditions 


Yield  Stress 

aT0 

CFT90 

ac0 

a  90 

c 

O  45 

CTxv 

rCo 

rc 

■  90 

Experimental  Values  (ksi) 

138.61 

130.9 

157.67 

135.44 

144.01 

74.76 

0.685 

0.507 

Normalized  Exp.  Values 

1.02 

0.97 

1.16 

1 

1.06 

0.55 

0.685 

0.507 

Predicted  Values 

1.04 

0.94 

1.14 

1.03 

1.04 

0.58 

0.69 

0.501 

0.125  in 


Figure  1 .  (a)  Specimen  geometry  for  the  quasi-static  experiments,  (b)  specimen  geometry 
for  the  experiments  using  the  Torsional  Kolsky  bar  set  up. 
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Figure  2.  Schematic  diagram  of  the  experimental  set  up  for  the  torsional  Kolsky  bar 

experiments. 
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Figure  3.  Schematic  diagram  of  the  channel  die  used  for  non-proportional  loading 

experiments 


Figure  4.  Yield  surface  generated  using  the  orthotropic  yield  criteria  from  Cazacu  et  al. 

(2005)  for  Ti-6A1-4V  titanium  alloy 


Symbols:  Experimental  data 
Line:  Orthotropic  model  correlation 
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Figure  5.  Correlations  using  the  anisotropic  yield  criterion  along  with  the  normalized 
experimental  yield  stresses  in  different  directions  in  the  plane  of  the  titanium  alloy  plate 


Figure  6.  Model  predictions  with  the  observed  shear  response  of  the  titanium  alloy 
converted  to  effective  stress  and  strain  during  monotonic,  strain-rate  jump  and  dynamic 
torsion  experiments  at  room  temperature  (296K).  Dynamic  and  0.989  s'1  effective  strain 

rate  predictions  are  adiabatic. 
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Figure  7.  Dynamic  shear  response  at  room  temperature,  followed  by  dynamic 
compression  response  at  room  temperature  of  the  titanium  alloy,  and  predictions  using 
KFIL  model  with  the  material  constants  determined  from  uniaxial  compression 
experimental  data. 


Figure  8.  Dynamic  shear  response  at  room  temperature,  followed  by  dynamic 
compression  response  at  422°K  of  the  titanium  alloy,  and  predictions  using  KFIL  model 
with  the  material  constants  determined  from  uniaxial  compression  experimental  data. 
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Figure  9.  Dynamic  shear  response  at  room  temperature,  followed  by  dynamic 
compression  response  at  588°K  of  the  titanium  alloy,  and  predictions  using  KHL  model 
with  the  material  constants  determined  from  uniaxial  compression  experimental  data. 


Figure  10.  True  stress-strains  for  the  non-proportional  loading  channel  die  experiment  at 
10'1  s'1  strain  rate  in  the  loading  direction.  Note  the  stress  starts  to  increase  in  the 
direction  2  (constrained  direction)  from  zero  once  the  specimen  comes  in  contact  with  the 
die  walls  at  approximately  12%  strain  in  the  loading  direction. 


Figure  11.  Quasi  static  response  of  the  titanium  alloy  under  non  proportional  biaxial 
loading  at  room  temperature  converted  to  effective  stress  and  strain  based  on  the 
orthotropic  criteria  and  predictions  using  KHL  model  with  the  material  constants 
determined  from  uniaxial  compression  experimental  data. 
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